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ABSTRACT 

Aims. To investigate the ingredients, which allow star-forming galaxies to present Lyce line in emission, we studied the kinematics 
and gas phase metallicity of the interstellar medium. 

Methods. We used multi-object near-infrared spectroscopy with MagellanfMMIRS to study nebular emission from z — 2 - 3 star- 
forming galaxies discovered in three MUSYC fields. 

Results. We detected emission lines from four active galactic nuclei and 13 high-redshift star-forming galaxies, including Ho- lines 
down to a flux of (4± 1)E-17 erg sec -1 cirr 2 . This yielded seven new redshifts. The most common emission line detected is [OIII]5007, 
which is sensitive to metallicity. We were able to measure metallicity (Z) for two galaxies and to set upper (lower) limits for another 
two (two). The metallicity values are consistent with 0.3 < Z/Z Q < 1.2, 12 + log(0/H) ~ 8.2-8.8. Comparing the Lyce central 
wavelength with the systemic redshift, we find Av Lto _[o///]5oo7 = 70-270 km sec -1 . 

Conclusions. High-redshift star-forming galaxies, Lyor emitting (LAE) galaxies, and Ho- emitters appear to be located in the low mass, 
high star-formation rate (SFR) region of the SFR versus stellar mass diagram, confirming that they are experiencing burst episodes of 
star formation, which are building up their stellar mass. Their metallicities are consistent with the relation found for z < 2.2 galaxies 
in the Z versus stellar mass plane. The measured AvLy a -[oiii]5007 values imply that outflows of material, driven by star formation, could 
be present in the z ~ 2 - 3 LAEs of our sample. Comparing with the literature, we note that galaxies with lower metallicity than ours 
are also characterized by similar Avz, va -[0///]5oo7 velocity offsets. 

Strong F([OIII]5007) is detected in many Lyce emitters. Therefore, we propose the F(Lya)/F([OIII]5007) flux ratio as a tool for the 
study of high-redshift galaxies; while influenced by metallicity, ionization, and Lya radiative transfer in the ISM, it may be possible 
to calibrate this ratio to primarily trace one of these effects. 

Key words. Techniques: spectroscopy - Galaxies: star-forming - Galaxies: Lyman Alpha Emitters, Active Galactic Nuclei 



1. Introduction 

Characterizing star formation at high redshift is critical for the 
investigation and understanding of galaxy formation and evo- 
lution. This task is often difficult, however, since distant star- 
forming galaxies can be quite faint and many of their character- 
istic diagnostic features redshifted into different observational 
regimes. The identification of high-redshift objects during spe- 
cific evolutionary stages amidst foreground sources generally 
falls into two categories, either photometric methods which sam- 
ple portions of the integrated spectral energy distribution (SED) 
or narrow-band imaging which select specific emission lines at 
specific redshifts. 

Photometric methods (e.g., the Lyman break galaxy, or LBG, 
technique, Steidel et al. 1996; the Lyman alpha decrement, or 
BX technique, Steidel et al. 2004) have been extensively used 
to identify large samples of high-redshift galaxies, and even 
characterize their physical properties via comparisons to theo- 
retical SEDs, generated by various initial mass functions, star- 
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formation histories, stellar masses, dust reddening and inter- 
stellar medium (ISM) metallicities (e.g., Shapley et al., 2003; 
Pentericci et al., 2007). The advantage of the photometric tech- 
niques is that faint sources, occupying wide fields-of-view and 
wide redshift slices, can be detected, while the disadvantage is 
the lack of redshift confirmation. 

A more efficient way to detect galaxies, at least in a narrow 
redshift range, is by using narrow -band filters covering a specific 
line wavelength at a specific redshift, such as Lya (e.g., Cowie 
& Hu, 1998). In principle Lya could be observed up to z = 10 
from the ground. Over the past decade, astronomers have built 
large samples of Lya emitters (LAEs) at z > 2 (e.g., PJioads 
et al., 2000; Gronwall et al., 2007; Ouchi et al., 2008; Nilsson 
et al., 2009). Narrow-band techniques have been used to detect 
Ha, [OII]3727 and [OHI]5007 line emitters at various redshifts 
(e.g., Thompson et al., 1994; Hayes et al., 2010b; Hatch et al., 
2011; Shim et al., 2011). Such surveys, however, often require 
spectroscopic follow-up to weed out lower redshift emission- 
line galaxies which contaminate the samples. New multi-narrow- 
band surveys, with narrow -band filters covering Lya, [Oil] 3727 
and Ha+[NII] at z = 2.2 (e.g., Nakajima et al., 2012a; Lee et al., 
2012), now have the advantage to reduce almost completely the 
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contamination from low-redshift objects and also to calculate 
emission-line ratios, sensitive to galaxy metallicity. 

In the context of the MUSYC survey (Gawiser et al., 2006), 
1 we amassed two samples of LAEs at z — 3.1 and z — 2.1. 
Their physical properties were studied by Gronwall et al. (2007), 
Gawiser et al. (2007), Lai et al. (2008), Guaita et al. (2010), 
Guaita et al. (2011), Acquaviva et al. (2012), and Bond et al. 
(2012), which found LAEs to be in the lower mass range of the 
mass distribution of star-forming galaxies, observed during ac- 
tive phases of star formation, and characterized by low dust con- 
tent and compact morphology. These are the ideal samples to 
constrain any evolution (or lack thereof) between redshifts of 
z * 3.1 and z * 2.1. Ciardullo et al. (2012) found z 3.1 LAEs 
with luminosities brighter than L(Lya) > 10 424 erg sec -1 had al- 
most twice the number density of z — 2.1 LAEs in the Extended 
Chandra Deep Field South. Also, the number of galaxies with 
EW(Lya) > 100 A at z — 3.1 was almost twice the number at 
z — 2.1. While z ~ 3 LAEs were observed to be consistent with 
zero dust content (Gawiser et al., 2007; Nilsson et al., 2009; Ono 
et al., 2010), the SED fitting of LAEs at z ~ 2 showed a mod- 
erate amount of dust (Nilsson et al., 2011; Hayes et al., 2010a; 
Guaita et al., 2011). 

The SED analyzes also showed that z — 3.1 LAEs are low 
mass (10 8 - 10 9 M ), young galaxies (20 Myr) in active phases 
of star formation (Gawiser et al., 2007; Lai et al., 2008), while 
z - 2.1 LAEs are instead characterized by more variety. For 
instance, some z — 2.1 LAEs have 'ages since star formation be- 
gan' that are compatible with a few hundreds of Myr and the sub- 
sample of more massive galaxies (10 9 5 M Q ) are characterized by 
redder colors (e.g., Gu2011). This was interpreted as an indica- 
tion that there are a wider variety of ISM mechanisms, due to the 
spatial and velocity distribution of dust and neutral Hydrogen 
at lower redshift, which could be responsible for the escape of 
Lya photons (with delayed episodes of star formation among the 
proposals; Shimizu & Umemura, 2010). A recent scenario was 
described in Acquaviva et al. (2012), who used a Monte Carlo 
Markov chain code (GalMC, Acquaviva et al., 2011) to derive 
physical parameters from the SEDs of our z — 3.1 and z — 2.1 
LAE samples. They proposed that higher-redshift LAEs were on 
average more metal-poor, even if older, than lower redshift ones. 
The investigation of intrinsic properties and direct metalllicity 
estimations for individual galaxy could be very useful to confirm 
this hypothesis. Our current work here looks to address two main 
aspects of this scheme via estimates of the gas phase metallicity 
and an investigation of ISM kinematics. Star-formation activity 
and supernova explosions could be responsible for pushing ma- 
terial out of the galaxies in the form of outflows, which conse- 
quently leave imprints in the profile of Lya emission ( Verhamme 
etal.,2006). 

The estimation of metal line abundances is important to in- 
crease the statistics of metallicity estimations in star-forming 
galaxies (SFGs) at high redshift. While Balmer emission lines 
can be modeled by stellar synthesis codes (Schaerer & de Barros, 
2011), a deeper investigation of the [OIII] emission relative to 
the Hydrogen lines is useful to allow better SED modeling of 
optical emission lines. Notably, strong [ODTJ5007 emission can 
comprise a significant fraction of the K band magnitude at z ~ 3 
(McLinden et al., 2011), making its characterization critical for 
the correct interpretation of broad-band photometric properties. 
Moreover, McLinden et al. (2011) and Finkelstein et al. (2011) 
recently showed that significant velocities offsets between Lya 
and systemic redshift as measured by [OIII] can be present in 



LAEs, which are clearly important for an interpretation of the 
Lya emission. 

This pilot project focuses on obtaining rest-frame optical 
spectra for a sizeable sample of high-redshift galaxies via near- 
infra-red multi-object spectrographs (NIR-MOS). Only a few 
NIR-MOS instruments have been available to date. NIR spec- 
tra have been obtained by using the multi-object infra-red cam- 
era and spectrograph (MOIRCS) at the Subaru telescopes to 
observe cluster stars (e.g., Scholz et al., 2009), medium red- 
shift massive galaxies (e.g., Onodera et al., 2010; Hayashi et al., 
2011), BzK galaxies with K AB < 22 (Yoshikawa et al., 2010), 
Har emitters (Tanaka et al., 201 1), and quasars (Del Moro et al., 
2009). Recently the MMT and Magellan Infrared Spectrograph 
(MMIRS) was used to obtain NIR spectra of LAEs at z - 2.2 
(Hashimoto et al., 2012) with relatively bright optical and NIR 
magnitudes (at least compared to those observed here). 

We present here a new NIR spectroscopic survey, with five 
main aims: 

1) to demonstrate the capabilities for studying faint emission- 
line galaxies with MMIRS on the Magellan Clay 6.5m tele- 
scope; 

2) to confirm redshifts of photometrically-selected SFGs; 

3) to measure rest-frame optical emission line ratios and esti- 
mate metallicities of galaxies undergoing active phases of star 
formation; 

4) to place SFGs, previously detected via Lya or Ha emission 
lines, in an evolutionary sequence and investigate similarities 
and differences with respect to continuum-selected SFGs in the 
same redshift range; 

5) to probe the gas kinematics of LAEs in order to study chan- 
nels of the escape of Lya photons. 

We adopt AB magnitudes throughout, following the Two 
Micron All Sky Survey (2MASS) zero-point magnitudes, Kab = 

Kvega + 1-87. 

When we refer to the rest-frame wavelength of emission 
lines, we use vacuum-wavelengths of Ai ya = 1215.67 A, 
^[o///]5007 = 5008.239 A, A Ha = 6564.614 A. We denote the stel- 
lar reddening as E(B-V) and the nebular reddening as Eg(B-V). 
We assume a ACDM cosmology consistent with the WMAP 5- 
year results (Dunkley et al., 2009, their Table 2), adopting mean 
parameters of Q m = 0.26, Q A = 0.74, H = 70 km s -1 Mpc~', 
cr 8 = 0.8. 



2. Observations and sample definition 

MMIRS is a near-IR imager and multi-object spectrograph with 
a 4' x 6.9' field-of-view, 2 operating on the Magellan Clay tele- 
scope since the end of 2010. Mask design followed the MMT 
mask preparation procedure. 3 The objects studied in this paper 
come from two two-night observing runs in late 2010 and mid 
2011 (see Table 1). Weather conditions and instrument perfor- 
mance were more favourable in the first run. Two entire nights 
were allocated in Nov, while only half nights in May, 19th and 
20th. The half nights allocated to the second run resulted in air- 
masses that were generally higher than the first run, reaching up 
to airmass 2.0 on May 20th. 

We obtained spectra in three fields belonging to the MUSYC 
survey (Gawiser et al., 2006, hereafter Ga2006): two masks in 
the ECDF-S (Extended Chandra Deep Field South), one mask 



http://physics.rutgers.edu/~gawiser/MUSYC 



2 http://www.cfa.harvard.edu/mmti/rnmirs/Calibration/waverange_HK_l.png 

3 http://hopper.si.edu/wiki/mmti/MMTI/MMIRS/ObsManual/ 
MMIRS +Mask+Making 
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in the EHDF-S (Extended Hubble Deep Field South), and one 
mask in the SDSS-1030 region. 

In Table 2 we show the number of detected sources versus 
the targeted ones. In Table 3 we list the 17 sources from which 
we detect emission lines. The names designed in column 1 will 
be used hereafter. When filling the masks we chose three main 
categories of sources: 

1) Narrow-band selected LAEs from Gronwall et al. (2007; 
hereafter Gr2007) at z — 3.1 and in Guaita et al (2010; hereafter 
Gu2010) at z — 2.1. Also, we considered Lya and/or Ha emitters 
detected by Hayes et al. (2010) in their double blind (DB) sur- 
vey. The detected sources originally belonging to this category 
and studied here are LAE27, z3LAE2, z2LAE2, z2LAE3, DB8, 
DB12, and DB22 (Table 3 notation). 

2) SFGs as defined in Ga2006 and Guaita et al. (201 1), and 
studied in Berry et al. (2012), selected using the broad-band 
color drop-out technique at z ~ 3 (LBGs) and Lya- decrement at 
z ~ 2 ('BXs'). The detected sources originally belonging to this 
category and studied here are BX1, BX4, BX10, BX14, BX16, 
and LBG11 (Table 3 notation). Our spectroscopy revealed the 
presence of Lya emission with equivalent widths larger than 
20 A for BX1 and BX10, confirming them as LAEs. Also in- 
cluded in this category were BzK galaxies either from Blanc 
et al. (2008) with Kab < 22 or, in the case of Bzk25, se- 
lected from deeper VLT infrared spectrometer and array camera 
(ISAAC) imaging with fainter continua (Mark Dickinson and 
Jeyhan Kartaltepe, private communication). The latter was not 
expected a priori to present any continuum in the NIR. 

3) Active galactic nuclei (AGNs) expected to lie in the range 
2 < z < 3 and show strong rest-frame optical lines, selected from 
optical (Francke et al., 2008), radio (Miller et al., 2008), or X- 
ray catalog (Xue et al., 2011). The AGNs originally belonging 
to this category and studied here are AGN5 and AGN26. 

We designed the MMIRS masks giving first priority to 
Lya/Ha emitting galaxies observed at z - 2 and z — 3. 

Galaxies selected to be either Lya and/or Ha emitters were 
expected have faint continua and only display a few emission 
lines in the NIR spectra like [OIII] and Ha due to star forma- 
tion. NIR spectra of DB sources, where Ha fluxes were mea- 
sured photometrically, are very important for setting upper limits 
on Ha detection. We filled in the remaining space on the masks 
using SFGs and AGNs with or without firmly known redshifts 
but suspected to lie at z > 1.5. We also included low-redshift 
galaxies, which were expected to yield strong continuum and 
thus be useful tests of MMIRS' capability, as well as slit stars 
for mask alignment phase and flux calibration. In each of the 
four masks we included 3-4 square boxes with alignment stars, 
although during the second run we used three additional stars for 
flux calibration and additional alignment checks. 

Our z - 2. 1 and z - 3. 1 LAE samples had a median Lya flux 
of 4.0 x 10~ 17 erg sec" 1 cirT 2 (Gr2007, Gu2010). The Lya/Ha 
ratio can be significantly reduced from the theoretical value 
of 8.7 for case B recombination (Osterbrock, 1989) even by a 
small amount of dust reddening. To produce a reasonable esti- 
mate of the reddening expected, we assumed a Calzetti law with 
E(B-V)=0.04, which is typical of our sample LAEs (Ga2007, 
Gu2011), and that Lya photons experience the enhanced dust 
attenuation of nebular gas but no further amplification of dust 
column from radiative transfer. This predicts a Lya/Ha ratio of 
~3, implying a median Ha flux for our sample of 1.3 x 10" 17 
erg sec -1 ctrT 2 . Scaling from the predicted S/N per resolution 
element for the current HK grism (MMIRS S/N calculator, 4 ) 
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we anticipated a S/N of 6 in 4 hours for the lines predicted in 
our typical LAEs. We obtained total exposure times of 4.8 and 
4.7 hours for our two ECDF-S masks, and 3.2 and 2.7 hours for 
our SDSS and EHDF-S masks, respectively. Spectra were taken 
through the HK grism and H + K filter covering 1.25-2.45 ymx 
wavelength range, through a slit with width equal to 0.7" (0.5" 
for a few sources in ECDF-S masks). The filter transmission has 
the same throughput of 70% at 16000 and 21000 A. At 19000 A 
the throughput has its maximum, but the atmospheric absorption 
is maximum at this wavelength as well. 

During the Nov run we chose a dithering pattern of 0, 3, -2, 

2, -3 arcsec along the slit length to perform a better sky subtrac- 
tion. However, the shift of pair subtracted images could cause 
an overlap of close-by slits (Sect. 3). Also, a 3" dither caused 
alignment stars to move outside of their boxes in a few cases. 
Therefore, we adopted a 0, 2, -2 pattern in the May run. We re- 
aligned the masks frequently during the nights, typically after 
a series of 5-10 science frames. We observed telluric standard 
stars 5 2-3 times per mask per night, placing the telluric in 5 slits 
across the mask to achieve full spectroscopic coverage. 

MMIRS is characterized by a 0.2012"/pixel scale. It corre- 
sponds to a 3.5 (2.5) pixel resolution element for a 0.7" (0.5") 
slit width. For a sampling of 7 A/pixel, this means a line resolu- 
tion of 25 (18) A and a resolution R=800 (1100) at A = 20000 
A. 

3. Reduction 

A full account of our reduction procedure is reported in the 
Appendix A, while here we describe only the main details. 

We used the mmfixen code provided by the MMIRS instru- 
ment scientific team 6 to collapse all the information enclosed 
in the multi-extension files. The output became a new multi- 
extension file, in which the 1st extension was the science frame 
normalized to counts/sec. 

To calculate the wavelength calibration map we used the 
COSMOS 7 public software package, which was originally de- 
signed to reduce optical multi-object spectra obtained with the 
Magellan IMACS and LDSS instruments and updated to work 
with NIR MMIRS data. Positional files (*.SMF) were generated 
for each mask by the MMIRS team (and mmirs2smf.py code) to 
make MMIRS mask frames compatible with COSMOS. We ob- 
tained a rough wavelength map using COSMOS align-mask and 
ad just-offset codes and then a final solution running COSMOS 
ad just-map. 

The most prominent features in our NIR spectra were the at- 
mospheric sky lines. They were used as a reference for the wave- 
length calibration, but afterwards they were eliminated from the 
science frames using an ABBA source offset procedure. We first 
subtracted two frames, A and B, adjacent in time to remove 
the bulk of the constant and variable IR backgrounds, obtain- 
ing background-subtracted (A-B) or (B-A) frames. Depending 
on the sky variation at a specific moment of the night and air- 
mass, in the A-B and B-A frames sky lines are reduced by about 
50-80%. We then applied offsets to align the dispersed source 
spectra from the A and B frames and summed the signals of 
aligned frames of the same mask for the entire run. The shift was 

5 http://hopper.si.edu/wiki/mmti/MMTI/MMIRS/MMIRS+Pipeline/ 
Telluric+Correction 

6 http://hopper.si.edu/wiki/mmti/MMTI/MMIRS/ 
MMIRS +Pipeline 

7 http://hopper.si.edu/wiki/mmti/MMTI/MMIRS/ 
MMIRS +Pipeline/MMIRS +COSMOS 
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determined based on the input dithering pattern. The remaining 
sky line residuals were due to rapid temporal sky line variations 
and variations along the slit for individual lines. These variations 
dominate the background noise. 

We then applied the previously calculated wavelength map 
to the completely reduced frames. This was performed using the 
routine COSMOS extract-Id spec, which applies the wavelength 
map and and splits a multi-slit frame into individual slit 2D spec- 
tra. We used the IRAF twodspec.apextract.apall task to extract 
ID spectra. 

Giant AO stars (without strong features in their NIR spectra) 
were used as flux standards to generate the sensitivity function 
needed for flux-calibrating science frames, as well as to correct 
telluric absorption features. 

We chose extraction a pall apertures of +3 (+5) pixels for 
extracting ID standard and science source spectra, respectively. 
The fractional missing flux due to the wings of the PSF falling 
outside the extraction aperture (not to be confused with slit 
losses) is estimated to be less than 10% in both the runs. 

In the case of the Nov run, we compared our flux calibration 
for emission lines from AGN5 to those reported in the literature 
by van Dokkum et al. (2005); see Sect. 5 for more details. In the 
case of the May run, we compared the obtained sensitivity func- 
tions for our slit stars with their literature magnitudes, finding 
good agreement. 

4. Emission line fitting procedure 

After extracting ID spectra, we fitted emission lines and con- 
tinua with the python scipy optimized. leastsq function. The 
parameters of the fit were central wavelength (/lo), dispersion 
(o~ gauss ), and amplitude (A gauss ) of a Gaussian curve (f^ = 

(A gauss / V27r) exp[-(/i - A ) 2 /2a 2 auss ]) for emission lines and in- 
tercept of a flat continuum. 

To easily identify the first-guess fitting parameters we 
smoothed all raw spectra by a ~30 A box. We performed a lin- 
ear fit to the smoothed continua at ~ 16000 A and -21000 A, 
excluding emission line wavelength regions from the continuum 
fit. We fitted the smoothed-spectrum continuum-subtracted Ha 
and [OIII]5007 emission lines with single Gaussians to measure 
Aq (hence the source's redshift) and the first-guess <T gauss and 

A gauss . 

We performed the fit a second time to obtain best- 
fit o-gauss (o- gaU ss,best-fit) values for three Gaussians fitted to 
Ha+ [Nil] doublet simultaneously, fixing the redshift inferred 
from Aq in the previous step and varying cr gauss and A gauss around 
the first-guess values. We fixed the ratio of the [Nil] lines to 
1:3 (Osterbrock, 1989). Even if [Nil] is not clearly seen in 
the 2D spectrum, this procedure was important to estimate the 
contribution of [Nil] blended to Ha. Atek et al. (2011) found 
that the [NH]doublet contribution to Ha was on average 8% in 
their sample. We similarly fitted the [OIII]doublet + H/3 with 
three Gaussians, where we fixed the [OIII]doublet ratio to 1:2.98 
(Storey & Zeippen, 2000). We excluded pixels strongly contam- 
inated by sky-line residuals from the fit. 

After obtaining the best-fit Ao and cr gauss ^est-fit for the six 
emission lines in the smoothed spectrum, we fixed these val- 
ues and searched for the best-fit A gauss (A ?fl „ JS /, e „_j/ f ) in the raw 
spectrum. From this last amplitude estimation we calculated the 
integrated line flux (F(\me)=A gausSybesr - fi , x o- gmlssM ^ fit x V27r). 
We estimated the lcr uncertainty on the Aq, Cgaussjest-fiu an d 
A gaU ss,best-fit values assuming a single parameter of interest, cal- 
culating the diagonal elements of the fit covariance matrix. 



The [OIII]doublet was a robust feature to find because it is 
composed by two separate emission lines with a known flux 
ratio. If we saw [OIII], but not Ha, we assumed an Ha upper 
limit of twice the continuum root-mean-square (rms) error at the 
wavelength of Ha multiplied by the MMIRS resolution of 25 A. 
The same method was applied in the case we detected Ha but 
not [OIII]. 



4.1. Uncertainties of the emission line parameters 

To estimate line integrated flux uncertainty, we used the parame- 
ter uncertainty derived from the best fit covariance matrix (Sect. 
4). However, additional factors could increase the total error bud- 
get. 

We performed a few tests to verify the goodness of the re- 
duction and calibration processes. We calculated the observation 
point spread function (PSF) by measuring the median full width 
half maximum (FWHM) of the telluric standard spectral pro- 
file. Because telluric standards were observed in 5 slits along the 
mask spatial axis, we could estimate the average PSF variation 
along the mask. In addition to this, we could estimate the pos- 
sible variations of the FWHM spectral profile with time in the 
series of 5 sub-sequent exposures. This could be related to the 
alignment, which might have varied with time. But in a series of 
five 1 sec-long exposures, AFWHM was just 0.1". We did the 
same test considering the spectral profiles of stars in the science 
masks as well. The FWHM and the AFWHM are comparable to 
the ones measured for telluric stars. In the Nov run, for example, 
FWHM=0.8"±0.2" for telluric and science frames. This means 
that our flux calibration could compensate for slit losses, in the 
case extraction aperture and slit width were the same for telluric 
and science spectra. The fraction of the flux in the wings of the 
spectral profile ending outside the extraction aperture should be 
less than 10%. 

Flux calibration uncertainties are estimated to be less than 
10% in the Nov masks and up to 30% in the May masks. This 
was produced by the standard deviation of the flux among stan- 
dard spectra taken from the same slit along the run, due to 
changes in weather conditions. 

Sky line residuals could also affect the integrated flux estima- 
tion. We avoided significant pixels with residual emission from 
the sky in the line fitting. In the case of [OIII], the presence of 
a doublet and a fixed ratio 2.98:1 compensated for the potential 
increase/decrease of flux due to sky line residuals. 

Depending on the night seeing/slit width, sky emission lines 
can be broadened, producing an uncertainty in the wavelength 
calibration and so in the observed line center as well. The wave- 
length map needed to wavelength calibrate our spectra was suc- 
cessfully generated through the COSMOS package, but an un- 
certainty of 10 A is the upper limit in the wavelength uncertainty 
due to broadened sky lines. To estimate this value, we extracted 
sky spectra from random slits of a few science frames, in which 
the ABBA procedure was not applied yet. We also average com- 
bined all these science frames the same way the final reduced 
ones were combined and extracted the sky spectrum from the 
combination. Then we compared the emission line wavelengths 
with the tabulated ones. 

All these uncertainties are included in central wavelength 
and integrated line flux total error budgets, reported in the fol- 
lowing figures. 
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5. Sources detected in MMIRS spectra 

In Table 3 we list the 17 sources, for which we detected emis- 
sion lines (Ha and [OIII]5007 have signal-to-noise ratio, S/N 
> 3). AGNs, LAEs, and DB galaxies had the highest fraction of 
successful detections. The main reason for the low fraction was 
that the sources were generally faint in K band and that in some 
cases a particular emission line was overlapped with a strong sky 
line residual. 

Several of the targetted AGNs exhibited strong, broad emis- 
sion lines clearly seen in their 2D spectra (Fig. 1). Therefore, 
these sources were good candidates to test the sensitivity and ca- 
pability of MMIRS, as described below. In Fig. 2 we show the 
AGN spectra in black and we over-plotted the best fit Gaussian 
curve in the [OIII] and Ho- wavelength region as red curves. The 
blue spectrum is the sky spectrum from Rousselot et al. (2000). 
Fig. 3 shows the line flux ratio diagram from Kewley et al. (2001) 
where SFGs appear to be separated from AGNs (upper left side 
in each panel). We used this diagram (originally proposed by 
Baldwin et al., 1981) to either identify and/or confirm AGNs 
among our sources. 

5. 1 . Active galactic nuclei 

AGN5 was previously studied in (van Dokkum et al., 2005, here- 
after vD05), using Gemini/GNIRS. They used the GNIRS short 
camera with 32 lines mm -1 and 0.68" wide, 6.2" long slit. Their 
total exposure time was ~5500 sec, with a seeing of ~ 0.7". We 
used a 0.5" slit in the H + K grism configuration for a total ex- 
posure time of 17280 sec (4.8 hours). The estimated total spec- 
troscopic throughput is ~ 30% in K and ~ 35% in H. This af- 
forded us the chance to cross-check the performance of MMIRS 
against high-quality Gemini/GNIRS data. We mentioned before 
that our flux calibration was able to compensate for slit losses. 
By definition it also accounted for telescope aperture and sys- 
tem efficiency. However, for this object we needed to take into 
account an additional slit loss due to a narrower slit of 0.5" dur- 
ing a night with 0.8" seeing. This implied that we needed to in- 
crease the integrated flux by about 20% to compare with van 
Dokkum's results. As described in the previous section, we fitted 
the continuum-subtracted emission lines with Gaussian curves to 
obtain the integrated flux and its uncertainty, after correcting for 
slit losses. Sky emission line residuals at 16128 A were excluded 
from the fit of the [OIII] doublet and we obtained F([OIII]5007) 
= 32.5 ± 6.7 E-17 erg sec" 1 cm 2 , F([OIII]4959) = 10.9 + 2.3 
E-17 erg sec 1 cm 2 . Faint sky lines lying in the middle of the 
Ha+[NH] complex at 21 1 16 and 21 156 A might still affect the 
results. We measured F(Ha) - 28.1 ± 7.9 E-17 erg sec -1 cm 2 . 
The continuum rms in the [OIII] and Ha wavelength ranges are 
not significantly different. 

Considering a PSF variation (AFWHM) of the order of 0.2" 
during the run and a flux calibration uncertainty of 10%, the 
MMIRS line flux measurements are consistent with the ones 
obtained using GNIRS by vD05 (F([OIII]4959)=(10.0+1.0)E- 
17 erg sec" 1 cm 2 , F([OIII]5007)=(24.0+2.4)E-17 erg sec" 1 cm 2 , 
and F(Ha)=(24.0±2.4)E-17 erg sec" 1 cm 2 ). 

The squares in the upper and lower panels of Fig. 3 represent 
our estimations of line ratios for AGN5 (filled) and previous van 
Dokkum's (open). Even if H/3 has very low S/N in our spectrum 
and so [OIII]/Hy6 ratio is consistent with an upper limit, AGN5 
appear to be located in the AGN region of the diagram. 

BzK25 was originally selected as a BzK galaxy. A public 
optical spectrum was available from the GOODS VIMOS survey 




-2 -1 



log([SII]doublet/Ha) 

Fig. 3: Line ratio diagrams as presented in Kewley et al. (2001) (solid 
lines). The solid lines divide AGNs from Starburst galaxies. Symbols 
indicate AGNs found in our survey as explained in the text, the filled 
(open) squares are our (van Dokkum et al., 2005) ratio estimations for 
AGN5, the black triangle (pentagon) represents Bzk25 (AGN26), the 
red star z2LAE3. SFGs located in the region typical of Starburst s are 
presented for comparison: DB12 as a magenta circle, BX16 as a cyan 
star, and the two LAEs from Finkelstein et al. (201 1) as orange dots. 



(GOODSXRb_002_l_q4_79_l.spec, Quality A optical redshift 
of 2.1449, Balestra et al., 2010). 8 

This showed CIII]1909 and a weak Hell in emission. This 
spectrum might indicate it was an AGN, although CIV and NV 
were not formally detected (we note there is a faint line at the 
location of CIV which appears to have a P-Cyg profile; such fea- 
tures can be observed for SFGs). The MMIRS spectrum did not 
show any continuum, but showed [OIII] and Ha emission lines 
(S/N([OIII])=7, S/N(Ha)=10). The filled triangle in Fig. 3 rep- 
resents this source and supports the idea it is an AGN. Looking 
at Fig. 2 we can see sky lines are not significantly affecting Ha 
or [OIII]5007. 

AGN26 is represented by the open pentagon in Fig 3. The 
line ratio diagram, we chose to show here, indicates it is an AGN. 
The [OIII] doublet is not affected by sky emission lines, but sky 
line residuals could affect the [SII] integrated flux. The red star 
corresponds to z2LAE3. Even if it was initially selected as LAE 
through narrow-band technique in EHDF-S, it appeared to be 
an AGN based on Fig. 3 line ratios. It is the brightest source of 
the sample, with L(Lya) was equal to 10 43 5 erg sec -1 (F(Lya) 
^ 1.0E-15 erg sec" 1 cm 2 and EW{Lya) res ,-f ramt , = 74 A) and it 
was the brightest source of the LAE sample in EHDF-S. 

For comparison, we show BX16 and DB12, which fall in 
the Starburst region of the diagram, as well as LAEs from 
Finkelstein et al. (2011). However, we note that the most mas- 
sive and most Lya luminous galaxies are generally found near 
the border of AGN and Starburst regions. 



8 http://www.eso.org/sci/activities/garching/projects/goods/ 
spectroscopy.html 
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Fig. 1 : 2D spectra of AGNs detected in the 4 masks, all characterized by Kab > 20. The spectra are aligned in wavelength. Poisson noise due to 
sky emission line residuals can be seen as black vertical stripes. 
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Fig. 2: Observed ID spectra of AGNs in grey and best fit in red. From the top AGN5, BzK25, AGN26, z2LAE3. Vertical dashed lines represent the 
wavelengths of sky line residuals. The blue curve is the sky spectrum by Rousselot et al. (2000), which indicates the relative intensity of the OH 
emission lines. It has been scaled according to the panel axis to be visible. [OIII]doublet (Ha) wavelength ranges are presented in the left (right) 
panels. 



5.2. Star-forming galaxies 



We adopted the same procedure described in Sect. 4 and 5 for 
fitting continuum and emission lines of all the sources listed 
in Table 3 and showed in Figs. 4 and 5. The SFGs detected 
in our MMIRS spectra showed an average Kab magnitude of 
22.3. We were not able to detect either emission lines nor con- 
tinua from BzK galaxies fainter than K AB = 23 or photometri- 
cally selected LAEs with L(Lya') fainter than lfj 4294 erg sec -1 . 
The most prominent emission line observed is [OIII]5007 (typ- 
ical S/N = 3-5), which is very intense in SFGs (Schaerer & 
de Barros, 201 1) and easy to identify thanks to its doublet. In 
Table 4 we report line fluxes, redshift, observed-frame equiva- 



lent width, and luminosities. This information was used to char- 
acterize the properties of SFGs in our sample. 



5.2.1. Physical properties 

In the past six years our sample of SFGs at 2 < z < 3 was 
studied via photometry and optical spectroscopy. In Gu2011 
we fitted their stacked SED and showed that stellar mass and 
reddening were the only parameters we were able to constrain 
with high significance. We found that the median stacked SEDs 
were able to reproduce the median stellar E(B-V) and mass 
values of a sample, we defined as the typical ones. We ac- 
counted for the spread in the sample population by using the 
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boot-strapping method to evaluate the uncertainties. It was also 
pointed out that stacking increased the signal-to-noise of ground- 
based photometric SEDs and revealed typical properties of our 
LAE samples. We needed LAE sub-samples (based on rest- 
frame UV magnitude, color, Lya equivalent width, and bright- 
ness at observed-frame 3.6 /mi separations) to identify unusual 
galaxy properties different from the typical ones. 

Among GOODS VIMOS public spectra (Popesso et al., 
2009; Balestra et al., 2010), we found that BX1 and BX10 were 
characterized by EW rest -f rame (Lya) > 20 A and BX16 was char- 
acterized by Lya in absorption. Thanks to MUSYC VIMOS 
and IMACS rest-frame UV spectra, LAE27 and z3LAE2 were 
confirmed to be LAEs, while LBG11 showed EW rest -f rame 
8 A (Fig. 4). In Table 6 we reported Ha/HB, Lya/Ha, and 
Lya/[OIII]5007 line ratios for all the SFGs of our survey. The 
integrated line fluxes were obtained by fitting MMIRS spectra 
shown in Fig. 4 and 5. 

For the objects with available rest-frame UV spectra, a good 
(stellar) reddening estimation could be obtained from the B slope 
of he rest-frame UV continuum. This is shown as E(B-V)yg in 
Table 6. Following Calzetti et al. (2000), B is related to E(B-V) 
through the following equations, 

A(1600) = 4.85+ 2.31 x 8, (1) 
A(1600) = 4.39 x Eg(B - V), E(B - V) = c x Eg(B - V), (2) 

where A(1600) corresponds to the galactic absorption in 
magnitude at A. rest -t rame - 1600 A, Eg(B-V) is the reddening of 
nebular gas, and E(B-V) that of stellar continuum. The constant 
c corresponds to a factor 0.44 + 0.03 in Calzetti et al. (2000). 
However, the Calzetti law was calibrated at low redshift and 
there are recent works showing that factor could be just c — 1 
at z ~ 2 (Erb et al., 2006). Hayes et al. (2010a) assumed c = 1 
as well, because they noticed that c = 0.44 significantly over- 
predicts the SFR(Ha) with respect to SFR(UV) for their sample. 
For our sample, SFR(Ha) and SFR(UV) are within a factor of 
2. From the SED fitting results in Gu2011, we also noted that 
0.44 < c < 1 seemed to be more appropriate to reproduce SFRs, 
but we could not estimate one precise value. It is beyond the 
scope of this paper to investigate if c is closer to 0.44 or 1, so 
we assumed the Calzetti law entirely. Of course to dust-correct 
emission line fluxes we used Eg(B-V). 

An alternative way to estimate dust reddening is the Balmer 
decrement (F(Ha)/F(H/J) Miller & Mathews, 1972). However, 
the low signal-to-noise HB line makes this an estimation which 
is less reliable than the previous ones. We reported the Balmer 
decrements and the implied stellar color excess, E(B-V) J s_ Q ., for 
4 sources in Table 6. 

We considered the stellar dust reddening from stacked SED 
fitting as the most robust measure to use for dust reddening 
correction, when a rest-frame UV fit could not be performed. 
Therefore, we referred to the sub-sample definitions in Lai et al. 
(2008) and Gu2011 to infer stellar mass and E(B-V) from their 
SED fits. The categories, presented in those works and consid- 
ered for the current analysis are IRAC -faint {f^^,,, < 0.3/iJ), 
IRAC-bright (f 3 , 6/jm > 0.3/J), UV-faint (R > 25.5), UV-bright 
(R < 25.5) LAEs, 'BX' and LBGs. The z ~ 2 LAEs detected 
in our MMIRS survey all had SEDs consistent with those of 
UV-bright LAEs (log(M/M )=9.1 [8.9-9.4], where square brack- 
ets correspond to the 68% confidence level range of stellar-mass 
values). The ones at z ~ 3 all belonged to the IRAC-bright sub- 
sample (M« ~ 10 9 - 10 10 Mq). The stacked sub-sample SED 
stellar mass was then rescaled to match the R band magnitude of 



individual sources. The SED for individual objects, performed 
in Hayes et al. (2010a), addressed stellar mass and reddening of 
DB sources. In Table 6 we reported the physical parameters from 
the SED fitting as E(B-W) SED and log(M*/M ). Their error bars 
accounts for the spread in the sample population (Gu201 1). 

To infer star formation rate (SFR) of SFGs and LAEs, we 
applied the Kennicutt (1998) equation. 

F corr {Ha) = F ohs (Ha) x lO - 4 *^-^^*,) (3) 
S FR corr (Ha) = 1.9E - 42 x F corr {Ha) x AjtD l {z) 2 (4) 

Only for three sources (BX1, BX10, BX16) we were able 
to estimate E(B-V) J a (and so Eg(B-V)g), we used to estimate 
SFR eorr (Ha). On the other hand Ha flux was an upper limit 
for the BX10 spectrum and it translated into an upper limit of 
SFR eo „.(Ha). In the case of BX4, BX7, BX14, and DB12, the 
SFR(Ha) was assumed to be a lower limit of the SFR coff (Ha) 
given the uncertain correction factor. On average SFR C0)T (Ha) 
was found to be within a factor of two of the value estimated 
from the rest-frame UV. 

As an example of the quantities introduced above, we present 
here the calculation relative to BX1. It was characterized by a 
quality A spectrum in GOODS survey. Fitting the spectroscopic 
rest-frame UV continuum with a power law, f A ~ A 13 , we ob- 
tained B best fu = -1.6 + 0.11 and so E(B-V) /3 = 0.12 ± 0.02, 
Eg(B-V) / 3 = 0.27 ± 0.04. The Balmer decrement is measured as 
F(Ha)/F(H/J) = 2.6±2.5 i.e. F(Ha)/F(Hy8) < 5.1, which implies 
Ep- a < 0.63. Its stellar reddening is then E(B-V) /; _ Q , < 0.54, 
roughly consistent with the SED estimation, but not very con- 
clusive. 

BX1 was originally selected to have R = 24.65 and it be- 
longed to the UV-bright LAE sub-sample (the R band magnitude 
of the stacked UV-bright sub-sample SED is 25.0). We, there- 
fore, assumed the stacked SED shape was representative of this 
object, but we rescaled by a factor 1.4 the stellar mass to ac- 
count for the individual object having a brighter R mag, yielding 
log(M/M o )=9.25[9.05-9.55]. 

Assuming E(B-V) /3 = 0.12 ± 0.02 and a Calzetti law for 
dust attenuation, we calculated F corr (Ha) = (9.4±2.1)E-17 erg 
sec -1 cm 2 and SFR eo „(Ha) = 28.5±8.2 M yr _1 . This value was 
roughly consistent with the SFR averaged over the last 100 Myr 
found in Gu2011 for UV-bright LAEs (< SFR >i o=10[6-70] 
M Q yr~'). Following the Kennicutt (1998) equation for the SFR 
implied by the ultraviolet continuum, BX1 is characterized by 
an SFR co „.(UV, R band) around 50 M Q /yr. 

5.2.2. Line velocity offset 

In the attempt to measure SFG properties we also studied ISM 
kinematics. [OIII]5007 can be used as a proxy of the galaxy 
systemic redshift (McLinden et al., 2011; Finkelstein et al., 
2011). Fitting Lya (from either GOODS or MUSYC spec- 
tra) and [OIII]5007 (from MMIRS spectra) emission lines with 
Gaussian curves we obtained best fits for their central wave- 
lengths. To compare GOODS with our redshifts, we first applied 
the heliocentric correction at the time of the observations. It was 
v=+22.2 km sec -1 at Cerro Paranal and v=+16.7 km sec -1 at 
Las Campanas for BX1, as an example. Then, the offset velocity 
is estimated with the following equation (Table 6), 

Av L ya-[O///]5007 = C X ( - l) (5) 

\ 1 + Z[OIII] ) 

where rest-frame wavelengths are A Em = 1215.67 A and 
hoimswi = 5008.239 A. 
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Fig. 4: Observed-frame ID spectra (black curve) of the sources showing either Lycc in emission or in absorption in rest-frame UV spectra. In 
red we show the best fit model spectra. From the top BX1, BX10, LAE27, z3LAE2, BX16, LBG11. The first, second and third panel columns 
represent Lya (from MUSYC and GOODS-VIMOS surveys), [OIII]doublet and Hor wavelength ranges respectively. Dashed vertical lines indicate 
the wavelength of sky emission lines, which could leave residuals in the science spectra. Horizontal dashed blue lines indicate the continuum rms 
in the wavelength range where we did not see an expected emission line. As LAE27, z3LAE2, and LBG1 1 are at z ~ 3, Ho- line would fall outside 
the MMIRS coverage. The rest-frame optical emission lines are unresolved. 



5.2.3. Gas phase metallicity 

In Table 5 we present the line ratios used to infer metallicity, 
parametrized by 12+log(0/H). We calculated [Oni]5007/H^, 
[OIII]5007/[NII]6584, and N2 ([NII]6584/Ha), which relate to 
metallicity through the Nagao et al. (2006, their Fig. 8 and 17) 
and Maiolino et al. (2008) calibrations. [OIII]5007/[NII]6584 
is sensitive to reddening, so it has been corrected using either 
E /3 (B-V) or E iS£0 (B-V) values (when the rest-frame UV spec- 
trum is not available) and the Calzetti law (Table 6). We re- 
ported [NII]6584/Ha for six sources, but only for three of them 
(BX16, BX4, DB12) the ratio is bigger than its error. It is ex- 
pected for SFGs with negligible or no AGN contribution that 
[Nil] is a weak emission line, hence it is difficult to use it to 
constrain metallicity. Therefore, in addition to N2, we also de- 
termined metallicity from two other emission line calibrations. 
The two rows in the third column of Table 6 represent the values 



inferred by the two branches of [OIII]5007/HjS vs 12+log(0/H) 
relation. For comparison with the literature we also adopted 
the Pettini & Pagel (2004) calibration of N2 (second row of 
7th column) and 03N2 ([OIII]/Hy8 / [NII]/Ha). In the last col- 
umn we show as 12+log(0/H) a </ opfe ,</ the value we will use in 
the discussion section. Due to low signal-to-noise from most 
[Nil] and H/3 emission lines, we were only able to constrain 
12+log(0/H) only for two sources (BX16 and DB12) and to set 
lower (upper) limits for two (two) others. In the case of LBG1 1, 
[Olll]5QQ7/H/3 = 2.8 ± 2.4 gave 12+log(0/H)< 9.0, which is not 
informative. In the case of z3LAE2, [OIH]5007/Hy8 = 1.7 ± 0.8 
was consistent with 12+log(0/H) of [7.0-7.2] and [8.5-8.7]. We 
considered the upper branch of the [OIII]5007/HyS calibration the 
most representative of this source metallicity, because EW(HyS) 
was measured to be at least 10 A (Hu et al., 2009). 

In the range of line ratios we measured for our sample, N2 
calibration from Pettini & Pagel (2004) and Nagao et al. (2006) 
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Fig. 5: Observed ID spectra (grey curves) of the SFGs detected in our survey, without spectroscopic rest-frame UV counterparts (DB8, DB12, 
DB22, z2LAE2, BX4, BX7, BX14). [OIII]doublet and Ho- regions are shown. If only one of these is detected, red vertical dashed lines indicate 
the location of Ho- ([OIII]5007) lines which would be implied by the detected-line redshift for DB8, DB22, z2LAE2 (BX4). Dashed grey vertical 
lines indicate the wavelength of sky emission lines, which could leave residuals in the science spectra. Horizontal dashed blue lines represent the 
continuum rms in the wavelength range where we did not see an expected emission line. 



agreed within the error bars. As an example, [NII]/Ha' = 0.2 for 
DB 12. This ratio corresponds to 12+log(0/H)=8.7 or 8.5 follow- 
ing Nagao et al. (2006) or Pettini & Pagel (2004) respectively. 
This difference of 0.2 dex was within the line ratio error bar. 

As shown by Kewley & Ellison (2008), metallicity calibra- 
tions based on different approaches (such as line ratios and direct 
measurement of electron temperature) could produce a system- 
atic deviation of 12+log(0/H) up to 0.7 dex. However, the analy- 
sis of (Kewley & Ellison, 2008, their Fig. 2) showed that N2 and 
03N2 have the lowest deviation from all the other possible cali- 
brations at 9 < log(M»/M ) < 10, which is the range of masses 
covered by the sources in our sample. 

We chose recent literature estimations of metallicity of SFGs 
and LAEs to compare with in a meaningful way. 

Atek et al. (2011) used a direct method, measur- 
ing [OIII]4363 (dependent on electron temperature) versus 
[OIII]5007 to infer metallicity of line emitter sources at 0.35 < 
z < 2.3. This method is successfully applied to derive low- 
metallicities (7.1 < 12+log(0/H) < 8.3), because at higher 
metallicities [OIII]4363 becomes too weak to be measured. 
They discovered a very metal-poor source at z — 0.7 with 
12+log(0/H) = 7.47 + 0.11 and log(M H ,/M )< 8. Here, we took 
this source as an example of a very metal-poor line emitter. 
We also compared with Xia et al. (2012), which used the R23 
(([OII]3727+[OIII]5007)/Hy8) method to estimate metallicity of 
SFGs with strong emission lines at z > 1. 

In addition to them, we considered sources at z > 2. The star- 
forming galaxy, named BX418, was found to have 12+log(0/H) 
= 7.9 ± 0.2 (Erb et al., 2010, direct method). With EW(Lyar) 
= 54.0 + 1.2 A, it is one of the most metal-poor LAEs known 
at high redshift. Finkelstein et al. (201 1) placed upper limits on 
metallicity for two LAEs at z = 2.3 and 2.5, using N2 and 03N2 
from Pettini & Pagel (2004). They can be directly compared to 



our sources. Nakajima et al. (2012b) selected a sample of LAEs 
at z ~ 2.2 and were able to set a lower limit of 12 + log(0/H) > 
7.9 on the stacked object, using a combination of [OII]/Hy6 and 
N2 calibrations. Nakajima et al. (2012a) used N2 from Maiolino 
et al. (2008) to infer 12+log(0/H) of individual sources, there- 
fore their results can also be directly compared with ours. 

We could apply neither R23 nor direct methods because in 
our MMIRS spectra we did not cover [OIII]4363 or [OII]3727. 

In Fig. 6 we show the mass-metallicity (M*-Z) relations pre- 
sented in (Kewley & Ellison, 2008, their Fig. 2), with their same 
color coding. The shaded regions represent the ranges in stellar 
mass and metallicity for the literature results described above. 
Kewley & Ellison (2008) pointed out that their M*-Z relation, 
corresponding to the direct method, was characterized by high 
statistical uncertainty. 

At the stellar mass of BX418 (log(M/M ) = 9.0), the dis- 
crepancy between the direct and N2 methods is less than 0.2 
dex. In the range of masses covered by the Xia et al. (2012) 
sample, it is 0.3 dex between N2 and R23 estimators. In the 
range of masses probed by Nakajima et al. (2012b), it is less 
than 0.4 dex among all the calibrations. These discrepancies 
were within the error bars of the metallicity estimations of our 
sources. Therefore, our results could be compared with the liter- 
ature ones just mentioned, reasonably without scaling to a com- 
mon calibration (Sect. 6). 



6. Discussion 

We found the [OIII]4959, 5007 emission line doublet to be the 
strongest feature in MMIRS spectra of z ~ 2 - 3 sources. We 
detected the [OIII] doublet in 12 out of 13 galaxies, while we 
detected Ha in just 6 out of 10 z ~ 2 galaxies. This may be 
surprizing at first glance. However, for SFGs characterized by 
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Fig. 6: Mass-metallicty relations as presented in Kewley & Ellison 
(2008) for log(M/M Q ) > 8.5. The color-coding is the same as in 
their Fig. 2. The cyan curve is obtained from the R23 calibration, 
blue, dashed black, and dashed red from a combination of R23 and 
[OIII]/[OII], solid black from a combination of line ratios, such as 
[Oil], 11/3, [OIII], Ho-, [Nil], [SII], solid green from 03N2, solid yellow 
from N2, and solid red from the measurement of the electron tempera- 
ture, using [OIII]4363, 4990, 5007. Shaded regions indicate the mass- 
metallicity area occupied by literature results within the error bars. Also 
shown is the Atek et al. (2011) source in brown at the lowest mass 
range, Finkelstein et al. (201 1) upper limits in orange, Nakajima et al. 
(2012b) lower limit in green, and the range of objects from Xia et al. 
(2012) z > 1 in grey. The sources in our survey are characterized by 
9 < log(M/M ) < 10, where the N2 and 03N2 agree with all the other 
calibrations within 0.4 dex. 



slightly sub-solar metallicity, [OIII]5007 is predicted and mea- 
sured to be very strong (Schaerer & de Barros, 2011). Kakazu 
et al. (2007) conducted a survey to search for ultra-strong emis- 
sion line galaxies at < z < 1. Their results showed that 
these galaxies contribute roughly 10% of the measured SFR den- 
sity at that epoch. All the Ha emitters in their sample were 
metal-poor (Z < 0.45 Z Q ) galaxies, characterized by strong 
[OIII]4363, while [OII]3727 was found to be very weak com- 
pared to [OIII]5007. The higher-metallicity objects had, instead, 
stronger [OII]3727/[OIII]5007. They concluded that they were 
looking at galaxies different from local metal-poor dwarf galax- 
ies with larger [OIII]5007/[OII]3727 flux ratio and more mass. 
This suggested an early chemical enrichment phase in galaxies 
which are still growing. 

Also, in our spectra the sky background continuum can be 3 
times higher at 22000 A than at 16000 A, and we measured an 
average fractional background error at 22000 A twice the one at 
16000 A. In Fig. 4 and 5 we presented the spectra of the SFGs 
in our sample, separated into [OIII]doublet and Ha wavelength 
ranges. 

6.1. Star-forming galaxy evolutionary stage 

Thanks to our spectroscopic results, we were able to place SFGs 
in the evolutionary sequence proposed by Noeske et al. (2007) 
(SFR eorr (Ha) versus stellar mass relation defined for log(M/M Q ) 
> 9.5 at z < 2). It was discovered that SFGs form a distinct 
sequence with a limited range of stellar masses. This 'main se- 
quence' moves as a whole to higher star formation rate as red- 



shift increases. Therefore, galaxies of a given mass tended to 
present higher SFR at higher redshift; they were much more ac- 
tive on average in the past. One reason is the larger abundance 
of gas, depleted with time. For galaxies in the main sequence, 
the gradual accretion (as opposed to episodic and bursty events 
such as mergers) is the dominant process, and major mergers are 
not the dominant trigger of their star-formation activity and evo- 
lution. The outliers of the mass-SFR relation are galaxies with 
larger star-formation rate with respect to ordinary SFGs. 

To completely fill this relation, we should be able to estimate 
SFR(Ha) for low mass (M, < 10 9 Mg) objects at z > 2. But the 
objects we detected in our spectroscopic survey tend to be bright 
in K band and therefore to be more massive than few x 10 9 M . 
Ga2007 had first showed that z - 3.1 LAEs are, typically, 10 9 
M galaxies in their active phases of star formation and one of 
the highest specific SFR objects at that redshift. This is in favor 
of the idea that LAEs were still building stellar mass. If dust is 
an indication of a more evolved galaxy population and it is one 
of the reasons why Lya photons cannot escape a galaxy ISM, 
we expect LAE to be characterized by lower mass and higher 
specific star formation rate. Nilsson et al. (2009) found that z ~ 2 
LAEs instead presented a more massive (log(M/Mo) ~ 9.5) sub- 
sample. 

In Fig. 7 SFR corr (Ha) is plotted versus stellar mass for SFGs 
at z ~ 2 - 3. Small symbols correspond to the literature data we 
chose for comparison, as described in the previous section. The 
lines (solid) with associated scatter (dashed) represent the cor- 
relation derived for z ~ 2 SFGs in the GOODS field by Daddi 
et al. (2007) and from GMASS (Galaxy Mass Assembly ultra- 
deep Spectroscopic Survey) by Talia et al. (2012). Our LAE and 
SFG data points tend to be characterized by higher SFR than 
typical z ~ 2 SFGs. Atek et al. (201 1) observed that metal-poor 
galaxies at z ~ 1 - 2 are located towards the upper left cor- 
ner of Daddi's correlation, indicating that they are experiencing 
strong metal-poor burst episodes of star formation. The metal- 
poor source from Erb et al. (2010), the sources from Nakajima 
et al. (2012b), and the lowest mass LAE from Finkelstein et al. 
(201 1) are also located in the same region. 

In Fig. 8 we show metallicity versus stellar mass for our and 
the literature sources we chose for comparison. 

The sources detected so far in our survey are confirmed to be 
active star-forming, characterized by metallicity of the order of 
0.3 < Z/Z G < 1.2 (12 + log(0/H) ~ 8.2-8.8). As explained in 
Sect. 5.2.3, it is meaningful to compare our results with the ones 
presented in the figure. The black curves are the mass-metallicity 
relations from Maiolino et al. (2008), which used spectroscopic 
data from AMAZE (Assessing the Mass-Abundance redshift[-Z] 
Evolution). They estimated stellar masses assuming a Salpeter 
IMF, like in our SED fitting, and calculated the best fit param- 
eters for the relation at z=0.07 (dotted-dashed), z=0.7 (dashed), 
z=2.2 (solid, calibrated following Erb et al. (2006) observa- 
tions), and z=3.5 (long dashed). They showed evolution in the 
mass-metallicity relation stronger at lower stellar masses. Our 
six star-forming galaxy metallicity estimations are roughly con- 
sistent with Maiolino et al. (2008) M.-Z relation at z = 1-2. 
The Finkelstein et al. (201 1) source with stellar mass bigger than 
10 10 Md results in a very metal-poor object. 

Our UV-bright LAEs do not seem to be the most metal-poor 
objects at this redshift, even if for a couple of them we just set 
upper limits. This could suggest an early chemical enrichment 
during the growth phase of individual z ~ 2 - 3 Lya emitting 
galaxies (Kakazu et al., 2007). The Erb et al. and Atek et al. very 
metal-poor objects present values even lower than the lower limit 
inferred by Nakajima et al. (2012b) for z ~ 2.2 emitters. 



10 



L. Guaita: NIR spectroscopy of star-forming galaxies 



1000 i 




M/M 

Fig. 7: SFR corr (Ha) vs stellar mass for LAEs and SFGs at z > 2. Red and magenta colors indicate hya and Ho- emitting galaxies, while blue 
indicate UV-continuum-selected SFGs without Lya in emission. The blue star with cyan error bars represents BX16, which shows Lya absorption. 
The SFRCHa) upper limit is due to an upper limit in Ha flux, and lower limits are presented when it is not possible to estimate E i8 (B-V) from the 
rest-frame UV continuum slope for individual sources. The solid green line is the M.-SFR relation by Talia et al. (2012) for SFGs at z -2. The 
black line is the same relation for z ~ 2 SFGs in the GOODS survey (Daddi et al. (2007)). Dashed green and black lines indicate the lcr range of 
those relations. Measured values from Nakajima et al. (2012b) (green square), Finkelstein et al. (201 1) (orange dots), Atek et al. (201 1) (brown 
triangle), and Erb et al. (2010) (violet triangle) are also shown. 



The existence of a correlation between stellar mass and 
metallicity reflects the fundamental role that galaxy mass plays 
in galactic chemical evolution. It can mean that low-mass galax- 
ies could be already experiencing a chemical enrichment phase. 
Also, galactic winds could be responsible in removing metals 
from galaxies (Tremonti et al., 2004; Rodrigues et al., 2012). 
This hypothesis could explain the metallicity lower than 7.9 ob- 
tained by Finkelstein et al. (201 1) for a rare metal-poor source at 
that stellar mass. Tremonti et al. (2004) also noted that in large 
galaxies with deep potential wells, star formation activity and 
supernova explosions are not effective in pushing material in the 
form of an outflow. Inflows of metal-poor gas could either dilute 
the gas, reducing its metallicity, or turn on new star-formation 
episodes which would again metal-enrich the ISM. In addition to 
this, Mannucci et al. (2010) showed that the mass-metallicity re- 
lation, observed in the local and high-z Universe, may be due to a 
more general relation between stellar mass, metallicity and SFR. 
Based on their results, z ~ 2.5 SFGs with log(M t /M ) = 9.7 are 
characterized by 12+log(0/H) of the order of 8.5-8.6, for star 
formation rates of about 30 M Q yr _1 . For the same stellar mass 
and metallicity, the LAEs considered here could present even 
higher SFR. This can be an indication that different mechanisms 
could dominate in this kind of galaxies. We investigate the out- 
flow hypothesis in the next sub-section and how it is related to 
the Lya luminosity. 



We are not yet able to derive conclusions about either M„-Z 
or M*-SFR relations of LAEs and SFGs at the same redshift. On- 
coming spectroscopic data are needed to support the hypothesis 
that z ~ 2 LAE could be younger but more metal-rich than z ~ 3 
LAEs, as found by Acquaviva et al. (2012). 

6.2. Lya emitter ISM properties 

From the MUSYC (FORS, VIMOS, and IMACS; see Berry 
et al., 2012) and GOODS-VIMOS spectroscopic surveys, we 
were able to measure Lya fluxes (Fig. 9) for 2 LAEs at z ~ 2, 1 
SF galaxy with Lya in absorption, 2 LAEs at z ~ 3, and 1 LBG 
GOODS ID spectra were generated through an optimal extrac- 
tion by following the slit profile measured in each slit. This min- 
imized the slit losses for individual objects. Also, the GOODS 
team restricted the observations to airmass equal to 1.1, to min- 
imize the loss of light due to atmospheric refractions (Popesso 
et al., 2009; Balestra et al., 2010). 

The LyafHa flux ratio is generally investigated to study the 
ISM nature (Hayes et al., 2010a; Finkelstein et al., 201 1). In fact 
the presence of dust (to which the Balmer decrement is also sen- 
sitive) can significantly reduce the escape of Lyo- versus Ha pho- 
tons. For BX1 (Sect. 5.2.2) we calculated an escape fraction of 
Lyor photons of about 30% and for BX10 a lower limit of 7% 
due to F(Lya)/F(Ha) > 0.6. Finkelstein et al. (2011) calculated 
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Fig. 8: Metallicity estimation versus stellar mass. The big magenta circle, red and blue stars correspond respectively to DB, LAEs and SFRs with 
EW res( _/ rome (Lyor) < 20 A in our sample. Thin green upper limits and the green square are the data points from Nakajima et al. (2012a)(2012b), 
small grey dots are Xia et al. (2012) values for z > 1 emitting galaxies. Black long-dashed, solid, dashed and dotted-dashed curves are the M.-Z 
relations from Maiolino et al. (2008) for the labeled redshifts. The horizontal solid line indicates solar metallicity (12+log(0/H)=8.7). 



a ratio close to 7. First of all, they were observing two LAEs 
brighter in Lya than our originally continuum-selected galaxies. 
Also, besides the big uncertainty and the low statistics, we could 
be observing individual LAEs with a higher dust amount than 
theirs. But they calculated E(B-V) = 0.20 + 0.15 and E(B-V) 
< 0.13 from the Balmer decrement for their two objects. We es- 
timated E(B-V) = 0.12 + 0.02 (E(B-V) < 0.54) for BX1 andE(B- 
V) = 0.12 + 0.02 for BX10 from fitting the UV /? slope (Balmer 
decrement), which are consistent with theirs. Therefore, in ad- 
dition to dust reddening, other ingredients of LAE interstellar 
medium need to be considered to completely understand their 
properties, such as gas phase metallicity and kinematics. 

For most emission line objects, the Lya/[OIII]5007 flux ra- 
tio was calculated with high significance (S/N > 3). This ratio 
is affected by gas phase metallicity, kinematics, and Lya radia- 
tive transfer in the LAE ISM, related to the neutral Hydrogen 
content. Because metallicity could be a key ingredient which 
regulates Lya photon escape, we note that Lya/[OIII]5007, in 
addition to the Lya/Ha or together with SED fitting, can be in- 
vestigated to improve the understanding of star-forming galaxy 
ISM (temperature, metal abundances, geometry, kinematics). 
We noted a tentative trend between Lya and [OIII]5007 emis- 
sion lines, which could be related to the M,-Z relation. High 
statistics are needed to disentangle Lya/[OIII]5007 dependence 
on metallicity and dust content. 

We noted that SFGs at z ~ 2 - 3 show F(Lya)/F([OIII]5007) 
ratio between and 1, with z ~ 2 LAEs showing a ratio closer to 



1 (Fig. 9). In addition to the flux uncertainty reported in Table 4, 
the uncertainties described in Sect. 4. 1 are included in the [OIII] 
flux error budget. We used the public version of GalMC code 
(Acquaviva et al. (2011)), which makes use of Anders&Fritse 
2003 line ratios, to predict Lya/[OIII]5007 ratios of typical SFGs 
at different values of galaxy age and metallicity (Z). Assuming 
an age of 10 s yr, for 0.02 < Z/Z Q < 2.0, Lya/[OIII]5007 can be 
between 5 and 8, assuming case B recombination and no dust. 
For E(B-V) = 0.1 it becomes 2.5 < Lya/[OIII]5007 < 4.0 and 
for E(B-V) = 0.3 we calculated Lya/[OIII]5007 < 1.0, assuming 
the Calzetti law. Also, in the absence of dust [OIII]5007/Hy6 ~ 4 
and [OIII]5007/Ha ~ 1.4. Therefore, a Lya/[OIH]5007 < 1.0 
could be explained either with E(B-V) ~ 0.3, consistent with the 
SED fitting results, or involving radiative transfer effects, which 
preferentially absorb Lye with respect to Ha and [OIII] photons. 
However, the referred Finkelstein's LAEs present stronger Lya 
emission with respect to [OIII] than ours and significantly lower 
metallicity. 

It was proposed (Kunth et al., 1998; Heckman et al., 2001; 
Dijkstra et al., 201 1, as examples) that galaxy outflows, driven 
by star formation activity, are able to red-shift the central wave- 
length of the Lya photon /1-distribution. This way they appear 
invisible to neutral Hydrogen and can escape. To investigate 
the dependence of Lya/[OIII] on kinematics we measured Lya- 
[OIII] velocity offsets (Sect. 5.2.2). 

Using stacked rest-frame UV spectra of a few LAEs in our 
sample, Berry et al. (2012) measured velocity offsets of about 
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600 km sec between Lya and low-ionization absorption lines 
(LIS). As low-ionization absorption lines are thought to be trac- 
ing the neutral gas, this offset is thought to be produced because 
out-flowing neutral gas allows only the red side of Lya line pho- 
tons to escape, making it appear redshifted with respect to the 
systemic velocity. Here we estimated expansion velocities of 
possible outflows calculating velocity offsets between the cen- 
tral wavelength of Lya with respect to that of metal lines at the 
systemic redshift. 

In Fig. 9 we present the F(Lya)/F( [OIII] 5007) ratio ver- 
sus Avz, VQ ,_[o///]5007- In the figure we do not show the value of 
-897.8 ± 257.0 km sec 1 , calculated for the SFG with Lya in 
absorption (BX16). BX10 is characterized by 749.9 ± 188.3 km 
sec 1 and the highest uncertainty in [OIII] central wavelength 
(18095.5 ± 11.6 A). In fact this emission line is on top of sky 
line residuals (Fig. 4), which affect the line integrated flux and 
central wavelength. However, we can not exclude it being a high- 
velocity-offset source. For the other three LAEs at z ~ 2-3 (star 
symbols), the offsets are 70-270 km see -1 . These values are con- 
sistent with previous estimations from McLinden et al. (201 1) 
and (Finkelstein et al., 2011, updated by Chonis et al. private 
communication). For LBG11, we estimated Av Ljo ,_[ 0m ] 5 oo7 = 
344.5 ± 107.0 km sec -1 , which is consistent with the value esti- 
mated by McLinden et al. (2011) for their strongest [OIII]5007 
emission line LAE (F([OIII]5007)=(35.48+ 1.15)E-17 erg sec -1 
cm" 2 ) at similar redshift. The velocities we calculated here are 
also consistent with expansion velocities, Vexp = 100-300 km 
sec -1 , measured by Verhamme et al. (2008) for a sample of 11 
LBGs. Instead, Steidel et al. (2010) estimated stronger veloc- 
ity offsets of 445 + 27 km sec -1 for a sample of Ha emitters at 
z — 2.3. Their range of Vexp within 1 cr is plotted as a rectangular 
region delimited by dotted lines. 

While the outflow velocity, Vexp, is best estimated from the 
velocity difference between LIS, which trace the out-flowing 
material, and systemic redshift emission lines, AvL ya -[oiii]50cn 
can give upper limits for commonly assumed geometries. In 
Berry et al. (2012) and Shapley et al. (2003), they estimated off- 
sets between LIS absoiption lines and Lya, which could be up 
to three x Vexp. It was shown in Verhamme et al. (2008) that 
for an expanding shell model Aviya-us can be 3 (2) times of the 
expansion velocity, depending on the HI column density, NHI 
> (<)10 20 /cm 2 . The situation can change for different geome- 
tries. In the figure the range of Vexp values from Shapley et al. 
(2003) is presented as the rectangle delimited by the red dashed 
lines, rescaled to be compared to our Av Lw _[ ///] 5 oo7. The range 
of Vexp values from Berry et al. (2012) at z = 2 and z — 3 is also 
rescaled and it corresponds to the rectangle delimited by solid 
lines. 

We fitted Lya line profiles with single symmetrical Gaussian 
curves. The asymmetry of Lya profile of LAEs at z — 3.1 was 
quantified by McLinden et al. (2011) as the cr gauss of the red 
peak divided by the cr gauss of the blue peak (a). One of their two 
sources showed significant asymmetry in the Lya line profile, 
a — 2.1 ± 0.2 (a = 1.0 + 0.1 for the other source). They used the 
Hectospec multi-fiber spectrograph at the 6.5m MMT observa- 
tory, with an instrument resolution of ~6 A. We performed a test 
to estimate the centroid shift due to the symmetric versus asym- 
metric profile fit. A good fit to the data can be obtained with a 
symmetric profile and a central wavelength shift of about 1 A 
in the rest frame. It implies a velocity offset of 60-80 km sec -1 
for LAEs at z ~ 3-2, which is within the error bars of our 
measurements. 



We do not see any significant difference between redshift 2 
and 3, but the statistics are still too low to derive strong con- 
clusions. However, combining Fig. 8 and 9 we can observe that 
lower dust, lower metallicity, bright Lya- flux sources, such as the 
two LAEs studied by Finkelstein et al. (201 1), tend to be charac- 
terized by velocity outflows similar to those of our sources. The 
rare metal-poor object by Finkelstein et al. (201 1) could host an 
episode of chemical enrichment of pristine gas, also character- 
ized by a low amount of dust. The sources in our survey could 
be galaxies in more evolved stages, but with hundred km sec 
outflows as channels for Lya photon escape. 

We are also interested in understanding how galaxy physi- 
cal properties can relate to outflow expansion velocities. In Fig. 
10 we plot AvL ya -[onr\5007 versus stellar mass for our survey of 
Z ~ 2 - 3 SFGs and recent results described in the literature. The 
squares correspond to Hashimoto et al. (2012), who fitted indi- 
vidual galaxy SEDs to obtain stellar mass and used the MMIRS 
and Keck/NIRSPEC instrument to measure AvL ya -[oin]5om, 
AvLya-Ha, and to infer Lya velocity offsets. A weak trend of ve- 
locity offset with stellar mass can be seen, but individual cases 
of very high expansion velocity and log(M*/M )~ 9.0 need to 
be investigated in more detail. This trend could be explained 
considering that low mass galaxies are still in the process of 
building up their mass and have a weak gravitational potential 
to support strong outflows (Heckman et al., 2001; Orsi et al., 
2012). In this scenario log(M*/M )= 10 SFGs, like LBGs, could 
show higher velocity outflow. LBG11, for example, shows one 
of the highest Vexp of our sample. However, massive galax- 
ies can be characterized by high HI column density (NHI ~ 
10 2l) /cm 2 ), which would completely remove the Lya line blue 
peak and would leave a Lya emission line profile either asym- 
metric or double-peaked with the red side stronger than the blue 
one (Verhamme et al., 2008; Orsi et al., 2012). This would pro- 
duce a big Av Lw _[ 0// /]5oo7, without the need of strong Vexp. 
Multiple scattering can produce multiple peaks in the wings of 
the double-peaked profile (Kulas et al., 2012). 

Additional investigation of Lya line profile needs to be car- 
ried out to disentangle these hypotheses. 

7. Summary and conclusions 

In this paper we presented new spectroscopy of SFGs at z ~ 2 - 3 
covered by MUS YC, probing nebular emission. The original aim 
was to detect [OIII] and Ha emission lines in SFGs to better 
understand their ISM properties (Table 5) and to relate them to 
Lya emission line strength and profile, where observed. Details 
on MMIRS performance are presented in the Appendix. In Fig. 
2, 4, 5 we showed the extracted and fitted ID spectra for all the 
sources in our survey. In Table 4, 5 and 6 we presented measured 
line fluxes, flux ratios and galaxy physical properties. We list 
here the main conclusions of our NIR survey. 

1) We developed a reduction pipeline, which makes use of 
COSMOS software. We found COSMOS plus an ABBA source 
offset procedure successful in reducing MMIRS data (Sect. 3 and 
details in the Appendix A). We tested MMIRS behavior using 
high signal-to-noise spectra of AGNs in our masks. 

2) Through MMIRS spectra, we could confirm seven 
redshifts obtained via rest-frame UV spectra and determine 
seven(three) new redshifts of SFGs (AGNs). 

3) We estimated optical emission line intensity and line ra- 
tios for 13 SFGs, 6 of which were previously selected as either 
Lya and/or Ha emitters. Two additional SFGs were discovered 
to be LAEs in GOODS public survey spectra. We found the 
[OIII]5007 emission line to be the most prominent one in our 
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z ~ 2 - 3 galaxy spectra, which are all characterized by R < 
25.5. 

4) We adopted either rest-frame UV f3 slope or E(B-W) SED 
for reddening correction. We took advantage of the SED fitting 
results from Gu2011, Lai et al. (2008) and Hayes et al. (2010a) 
to estimate the stellar mass of the galaxies in our sample. In 
the M t -SFR(Ha) plane (Fig. 7), the SFGs in our sample were 
located towards the upper left corner of the correlation found 
for z ~ 2 SFGs by Daddi et al. (2007) and Talia et al. (2012). 
At fixed stellar mass, they presented larger SFR than typical 
SFGs. This implies that they are experiencing active bursts of 
star formation, which are building up their mass. One possibility 
can be the star-formation episodes implied by mergers (Noeske 
et al., 2007). As observed by Atek et al. (201 1), low-metallicity 
sources occupy the same locus of the M < ,-SFR(Ha') plane. We 
estimated metallicity from Nagao (2006) and Maiolino (2008) 
[Om]5007/HA [OIII]5007/[NII]6584, and [NII]6584/Ha cali- 
bration. We chose a few metal-poor sources from the literature 
to compare our metallicity estimations and demonstrated it is 
meaningful to compare them. Due to weak H/3 and [Nil] emis- 
sion lines, we could estimate metallicity for two SFGs (one of 
which is selected to be Ha emitter) and set limits for another 
four. At 9 < log(M„/M ) < 10 our sources have metallicity con- 
sistent with 12 + log(0/H) ~ 8.2 - 8.8. The two LAEs studied 
by Finkelstein et al. (201 1) were significantly brighter in Lya (3- 
6 times in the integrated line flux, Fig. 9), showed significantly 



lower metallicity, and had average lower SFR con (Ho') than ours. 
In the M,-Z plane, our sources agreed with the relations found by 
Maiolino et al. (2008) for z < 2.2 galaxies. The sources detected 
so far in our survey are confirmed to be actively star-forming, 
characterized by metallicity of the order of 0.3 < Z/Z Q < 1.2. 
However, we are not yet able to derive strong conclusions about 
the connection between M„-Z and M„-SFR relations of LAEs 
and SFGs at the same redshift. 

5) We were able to study nebular emission from four 
LAEs, one SFG with Lyo- in absorption, and one LBG with 
EW ofes _ /rame (Lya') = 34 A (EW resf _ /rame (Lya;) ~ 2 A). For high- 
redshift galaxies with both rest-UV and rest-frame optical spec- 
troscopy, our observations indicate that the two emission lines 
with highest signal-to-noise will often be Lya and [OHI]5007. 
For many dim objects, these will be the only lines detected. 
Hence, we find empirical motivation to form and study the 
Lyc/tOIII] flux ratio and to attempt to calibrate it both empir- 
ically and theoretically. This ratio is far from ideal as a physical 
tracer, being sensitive to a combination of metallicity, ionization, 
and Lya radiative transfer in the ISM, and it will be challenging 
to model. But given the arrival of a generation of multi-object 
NIR spectrographs and the difficulty of detecting Ha for sig- 
nificant redshift ranges of interest, the Lya'/[OIII]5007 ratio is 
worthy of investigation. 

Consistent with previous estimations AvL Va -[oni]5iKn was 
positive, indicating the presence of outflows. We found that 
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Fig. 10: Velocity offset between Lya central wavelength and systemic redshift versus stellar mass. The color coding is as in Fig 9. Green squares 
are from Hashimoto et al. (2012) for the sample of z ~ 2.2 LAEs, also studied in Nakajima et al. (2012a). The stellar masses of our sources are 
estimates from stacked SED fitting as explained in Sect. 5.2.1. 



LAEs are characterized by Vexp ~ 70-270 km sec -1 , in agree- 
ment with Berry et al. (2012). However, the scattering of Lya 
photons due to neutral Hydrogen atom, could mask the exact 
value of the outflow velocity. It was proposed by Tremonti et al. 
(2004) that outflows could remove metals together with neutral 
Hydrogen. Therefore, high-stellar-mass galaxies could have ex- 
perienced outflow phenomena in the past and now show lower 
metallicity. The outflow phenomena could still be going on. 
However, our sources present just slightly sub-Z . In Fig. 10 we 
showed velocity offset versus stellar mass. There is a weak trend 
for which higher mass SFGs could be characterized by higher 
velocity outflow. However, massive SFGs tend to be character- 
ized by NHI ~ 10 20 /cm 2 , which could produce a strongly red- 
shifted red Lya peak without the need for strong Vexp. 

The MMIRS spectrograph can be successfully used to study 
nebular emission of bright SFGs at z ~ 2 - 3, but a larger sample 
of galaxies is needed to derive strong conclusions. 
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L. Guaita: NIR spectroscopy of star-forming galaxies 
Table 1: Log of Nov 2010 and May 2011 runs 





Mask 


r vv nivn ) 


1 1 iimcl in rn l lis s ) 


t_exposure(t_ramp) 


N 


Nov, 1-2 


ECDF-S 1 






ls(ls), ls(ls), 3s(ls), 300s(5s) 


10, 5, 5, 10 


Nov 1-2 


fcdf-s i 




lamp 


ls(ls) 


10 


Nov, 1-2 


ECDF-S 1 




flat 




10 


Nov 1-2 


Pff)P C 1 


8 


NIR telluric standard. 


ls(ls) 


20 


Nov 1-2 


EfTVB Q 1 
CLUr-o 1 


8 


science( 1.13) 




58 


Nov, 1-2 


ECDF-S 2 






ls(ls), 300(5s) 


10, 10 


Nov 1-2 


FCDF-S 9 

1 A lyl . > Z, 




lamp 


IS[ IS) 


10 


Nov 1-2 


FCDF-S 7 

LLUI O Z, 




flat 


ls(ls) 


10 


Nov 1-2 


ELUr o z, 


8 


NIR telluric standard 


ls(ls) 


25 


Nov 1-2 


prrjp c 9 

DV-J-T-O Z, 


8 


science( 1 . 1 5) 




56 


May 19-20 


SDSS 




dark 


ls(ls), 3s(ls), 300s(5s) 


10, 10, 10 


May 19-20 


SDSS 




lamp 


ls(ls) 


5 


May 19-20 


SDSS 




flat 


3s(ls) 


5 


May 19-20 


SDSS 


1.0 


NIR telluric standard 


ls(ls) 


10 


May 19-20 


SDSS 


1.0 


science(1.2) 


300s(5s) 


38 


May 19-20 


EHDF-S 




dark 


ls(ls), 3s(ls), 300s(5s) 


10, 10, 10 


May 19-20 


EHDF-S 




lamp 


ls(ls) 


5 


May 19-20 


EHDF-S 




flat 


3s(ls) 


5 


May 19-20 


EHDF-S 


1.0 


NIR telluric standard 


ls(ls) 


10 


May 19-20 


EHDF-S 


1.0 


science(1.3) 


300s(5s) 


32 



Notes. The observations were carried out using the MMIRS HK grism and the H + K filter over 4 nights in three fields and four MOS masks. 
Column 1: Observation date. Column 2: Mask name. For a few sources of ECDF-S 1 and ECDF-S 2 masks we used a narrower slit of 0.5"; a 0.7" 
wide slit was used for the majority of the sources. Column 3: FWHM estimated from the continua of standard stars. The seeing was 0.8"(1") in 
Nov(May) run. Column 4: type of calibration and science exposures performed for each mask (airmass is average). Column 5: Individual frame 
exposures; the ramp time is also specified for each type of frame in parentheses (Appendix A). Column 6: Number of each type of exposure. The 
total exposure time of the science frames was 4.8 (4.7) hours for ECDF-S maskl (mask2), 3.2 (2.7) hours for SDSS (EHDF-S) mask. 



Table 2: Sources detected in MMIRS survey 



objects 


#targeted 


#detected 


#emission line detected 


#continuum detected 


LAE 


27 


4 


4 


1 


DB 


13 


3 


3 





BX 


29 


6 


6 





LBG 


1 


1 


1 





BzK 


6 


1 


1 


1 


AGN 


11 


2 


2 


2 


low-z 


5 


5 





3 


slit stars 


3 


2 





2 



Notes. In this table we show the number of objects detected versus the ones targeted for each category of objects in the masks. The object name 
corresponds to the original category of the sources, as explained in §2. 'low-z' are galaxies, whose BzK colors predict them to be at z < 1.4. They 
appeared bright in the observed H + K band and showed continua in MMIRS spectra. Together with slit stars these are useful to test the sensitivity 
and throughput of MMIRS. 



Table 3: List of sources showing emission lines in each mask 



name 


a 


S 


Rab 


Hab 


K A b 


literature z 


ECDF-S 1 














AGN5" 


3:32:31.46 


-27:46:23.20 


23.32 


21.18 


20.90 


z = 2.221 


BzK25'' 


3:32:50.83 


-27:48:04.68 




23.49 


23.05 




BX1" 


3:32:27.10 


-27:46:43.25 


24.65 






z = 2.223 


BX10" 


3:32:33.78 


-27:48:14.40 


24.76 






Z = 2.618 


LAE27 1 ' 


3:32:52.68 


-27:48:09.41 


25.51 


18.00 


18.21 


z = 3.085 


ECDF-S 2 














BX16" 


3:32:43.63 


-27:43:47.71 


24.03 






z = 2.317 


Z3LAE2' 


3:32:26.94 


-27:41:28.26 


24.13 






z = 3.114 


DB8 1 


3:32:36.14 


-27:42:23.07 






24.79 




DB12' 


3:32:31.52 


-27:43:51.38 


24.51 








DB22' 


3:32:41.32 


-27:45:49.13 


26.43 


24.88 


25.58 




SDSS-1030 














AGN26 


10:30:30.10 


5:21:06.19 


21.49 


20.72 


20.35 




BX4 


10:30:10.25 


5:22:15.77 


25.23 




22.79 




BX7 


10:30:13.39 


5:22:37.85 


23.87 


22.44 


22.48 




BX14 


10:30:18.86 


5:21:40.87 


24.18 


22.95 


22.33 




Z2LAE2 


10:30:07.60 


5:21:39.70 


25.02 




22.67 




EHDF-S 














Z2LAE3 


22:33:43.67 


-60:39:44.56 


24.25 








LBG 11 


22:33:22.54 


-60:38:49.63 


24.55 






z = 3.205 



Notes. Column 1: name of the source. Column 2 and 3: right ascension and declination. Column 4, 5, and 6: source R, H, and K magnitudes. 
Column 7: redshift if a z spec was known for that source. "Lyman Break and BX galaxies are selected following selection criteria described in 
Ga2006 and Gu2010. 6 radio catalog object, observed through a 0.5" wide slit. C DB source from Hayes et al. (2010). d BzK source from Mark 
Dickinson and Jeyhan Kartaltepe (private communication). e z — 3.1 LAEs from Gr2007. The sources in SDSS-1030 and E-HDFS masks all 
belong to the MUSYC survey. The R AB , H AB , K AB magnitudes were obtained from the references above. The literature reshifts came either from 
the references above or from the public GOODS spectroscopic catalog (see text). Tab. 4 lists our results. 
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L. Guaita: NIR spectroscopy of star-forming galaxies 
Table 5: Galaxy metallicity 



ID 

(1) 




-» 12+lg(0/H) 
(3) 


LOWJ9UU7 


-> 12+lgCO/H) 
(5) 




->12+lg(0/H) 
(V) 


[<«»|5l«)/ , |A'//|M 


-> 12+lg(0/H) 
(9) 


\2+\°(OIW) adopKd 
(10) 


Hp 

(2) 


I/V//16584 

(4) 


(6) 


(8) 


BX1 


7.9 ±7.4 


[all] 


7.9 ± 13.9 


> 8.2 


0.4 ± 0.7 


[all] 
<8.9 


20.3 ± 40.7 


>8.1 


> 8.2 


BX10 


2.8 ± 2.7 


[all] 
















BX16 


1.9 ± 1.0 


8.6[8.4-8.8] 

<7.3 


3.3 ± 2.9 


8.5[8.4-8.8] 


0.16 ±0.14 


8.6[7.9-8.9] 
8.5[8.3-8.6] 


12± 12 


> 8.2 


8.5+!!'? 

-0.1 


BX4 










0.3 ± 0.2 


> 8.4 
8.6[8.4-8.7] 






> 8.4 


BX7 


> 1.8 


[7.2-8.6] 


8.2 ± 8.3 


>8.1 


0.15 ±0.15 


< 8.5 

< 8.5 


> 11 


<8.3 


< 8.5 


BX14 






25.8 ± 82.0 


[all] 


0.05 ±0.14 


<8.7 
< 8.5 






< 8.6 


LBG11 


2.8 ± 2.4 


<9.0 








- 






< 9.0 


LAE27 


2.2 ± 2.5 


<7.7 
>8.1 








- 








Z3LAE2 


1.7±0.8 


7.1[7.0-7.2] 
8.7[8.5-8.7] 














8.6^;! 


Z2LAE2 


3.0 ± 1.3 


[all] 
















DB8 


3.0 ±2.6 


[all] 
















DB12 


2.9 ± 2.3 


[all] 


3.1 ±2.3 


8.5[8.4-8.9] 


0.22 ±0.15 


8.7[8.3-9.0] 
8.5[8.2-8.7] 


13.4 ± 14.4 


> 8.4 




DB22 


> 1.7 


[7.1-8.7] 

















Notes. Column 1: names of the sources as presented in Table 3. Column 2, 4, and 6: measured [OIII]5007/H/J, [OIII]5007/[NII]6584, and 
[NII]6584/Ho' emission line ratios. Column 3, 5, and 7: metallicity, parametrized by 12+log(0/H) implied by the ratios reported in columns 2, 
4, and 6 respectively. We used Nagao et al. (2006) calibrations for [OIII]5007/H/3, dust-corrected [OIII]5007/[NII]6584, and [NII]6584/Ha. We 
adopted the calibrations from Pettini&Pagel (2004) to infer [NII]6584/Hff (second row of 7th column) and 03N2 ([OIII]/H/3 / [NII]/Ha). Column 
9: the value of 12+log(0/H) we adopted for our sources as a combination of the information coming from the more conclusive ratios. In the case of 
z3LAE2 we considered the higher branch of [OIiI]5007/H/? calibration the most representative of its metallicity, because EW(H/3) was measured 
to be at least 10 A (Hu et al., 2009). 



Table 6: Galaxy properties 



ID 

(1) 


Ha 
Hp 

(2) 


Lya 
~Ha 

(3) 


Lya 
LO///]5007 

(4) 


Av (km/sec) 
(5) 


E(B-VV 

(6) 


E(B-VV_ ff 
(7) 


E(B-V) S£fl 
(8) 


log(M,/M ) 
(9) 


12+log(0/H) 
(10) 


BX1 


2.6 ±2.5 


2.5 ± 0.6 


0.8 ±0.1 


125.1 ± 108.1 


0.12 ± 0.02 


< 0.54 


O -29+O.U6 


9 25 +0 i 


> 8.2 


BX10 




> 0.6 


0.49 ±0.17 


749.0 ±205.8 


0.12 ± 0.02 




n 09+0.06 
U - J -0.23 


q 9 1 +0 30 




BX16 


1.7 ±0.6 






-897.8 ± 272.3 


0.35 ± 0.02 


< 0.01 


10+ U1U 

-0.30 


q q+0.1 
9> 


8 --'-o.j 


BX4 














35+°- U2 

"• JJ -0.04 




> 8.4 


BX7 














35+ 002 


9 6+ 01 


< 8.5 


BX14 


4.0 ± 3.2 










0.3 ± 0.5 


35+°- u - 


04+0.1 


< 8.6 


LBG11 






0.37 ±0.14 


344.5 ± 107.0 


0.07 ±0.14 




> 0.1 


> 9.95 


< 9.0 


LAE27 






0.08 ± 0.02 


167.8 ± 105.3 






< 0.1 


7 ' 7J -0.17 




z3LAE2 






0.11 ±0.03 


221.8 ±90.0 






n a 9 +0.O6 

"■ JM -a.2i 


q q-5+0.13 


8.6ti 


z2LAE2 














09 + 26 

-009 


9 r 03 




DB8 














47- ,uw 

"■ H ' +0.26 


9 76+" 1 




DB12 


4.4 ± 3.3 










0.42 ± 0.20 


26+ u ib 


54+01 




DB22 














22+ 26 


9 26+" 1 





Notes. Column 1 : names of the sources as presented in Table 3. Column 2, 3, and 4: Ha/H/3, Lya/Ha, and Lyay , [OHI]5007 line flux ratios. Column 
5: velocity offsets calculated using equation (5). Column 6, 7, and 8: dust-reddening estimations, obtained from fitting the rest-frame UV slope 
(E(B-V) /3 ), the Balmer decrement (E(B-V)^_ ff ) in the 2nd column, and from SED fitting (E(B-V)sed)- If continuum is too low to estimate /? slope 
we do not report E(B-V) /3 . BX16 spectrum provides F(Ho')/F(HyS) = 1.7 ± 0.6, i.e. F(Ha)[F(HfS) < 2.9 at 2<r, which is consistent with no reddening 
(E(B-V) / j_„ < 0.01). Column 9: stellar masses in unit of solar metallicity (M»/M ), inferred from the SED fit. Column 10: metallicity values taken 
from the final column of Table 5. 
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L. Guaita: NIR spectroscopy of star-forming galaxies 
Table 4: Emission line fluxes for all the sources 



Line 


jobs-frame 
central 


Zspec 


Flux 


rms f A (continuum) 




L 




A 




erg sec -1 cirT 2 


erg sec -1 cm -2 A -1 


A 


erg sec -1 


AGN5 


H/3 


15679.3+11.5 




(2.2±3.0)E-17 


32.1E-19 


>7 


(0.5+1. l)E+42 


[<?///]4959 


15988.5+12.9 




(8.9±1.6)E-17 


32.1E-19 


>30 


(3.3±0.6)E+42 


[0///]5OO7 


16143.6±4.6 


2.2234 ± 0.0009 


(26.6±4.7)E-17 / 


32.1E-19 


>80 


(9.9+1. 8)E+42 


[A7/]6548 


21100.8±9.7 




(7.2±1.2)E-17 


22.3E-19 


>30 


(2.7±0.5)E+42 


Ha 


21153.1+6.9 


2.2223 + 0.0011 


(22.5±5.9)E-17 / 


22.3E-19 


>100 


(8.6±2.2)E+42 


[A7/]6583 


21217.2+11.1 




(21.7±3.2)E-17^ 


22.3E-19 


>95 


(8.3±1.2)E+42 


[5//]doublet 


21651.8+2.7 




(8.7±5.2)E-17 


22.3E-19 


>40 


(3.3±2.0)E+42 


BzK25 


Lya(R) 


3830.3±0.5 




(24.3±2.0)E-17 


24.0E-19 


100 


(8.5±0.7)E+42 


H/3 


15303.0±6.0 




(7.2±4.3)E-17 


20.0E-19 


>18 


(2.5±1.5)E+42 


[0///]4959 


15605. 1+5.2 




(12.7±1.9)E-17 


20.0E-19 


>32 


(4.4±0.7)E+42 


[O///]5007 


15756.2+1.8 


2.1461 + 0.0003 


(37.7±5.7)E-17 


20.0E-19 


>94 


(13.2±2.0)E+42 


[JV//]6548 


20603.5+?? 




(2.8±1.3)E-17 / 


11.6E-19 


>12 


(1.0±0.5)E+42 


Ha 


20650.8+4.2 


2.1458 + 0.0006 


(25.1±2.6)E-17 


11.6E-19 


>100 


(8.8±0.9)E+42 


[A7/]6583 


20713.6+78.9 




(8.3±4.0)E-17 


11.6E-19 


>36 


(2.9+1. 4)E+42 


AGN26 


H/3 


14669.3+14.7 




(22.7±8.5)E-17 


27.1E-19 


130 


(6.8±2.5)E+42 


[0///]4959 


14958.9+12.6 




(18.0±2.5)E-17 


27.1E-19 


100 


(5.3±0.7)E+42 


[O///]5007 


15103.7+5.2 


2.0158 + 0.0010 


(53.5±7.5)E-17 


27.1E-19 


320 


(16.0±2.2)E+42 


[iV//]6548 


19739.8+36.3 




(20.5±3.6)E-17 


24.0E-19 


85 


(6.1±1.0)E+42 


Ha 


19784.7+24.1 


2.0138 + 0.0037 


(44.7±4.3)E-17 


24.0E-19 


180 


(13.4±1.3)E+42 


[A7/]6583 


19845.0+46.4 




(48.6±4.5)E-17 


24.0E-19 


200 


(14.5±1.3)E+42 


[5//]6718 


20252.0+11.6 




(5.3±3.6)E-17' 


24.0E-19 


20 


(1.6±l.l)E+42 


[5//]6730 


20288.1+17.5 




(16.1±7.3)E-17' 


24.0E-19 


65 


(4.8±2.2)E+42 


z2LAE3 


H/3 


14967.6+12.4 




(10.2±6.7)E-17 


72.0E-19 


>7 


(3.3±2.2)E+42 


[<?///]4959 


15263.1+6.5 




(44.7±2.7)E-17 / 


72.0E-19 


>30 


(14.4±0.9)E+42 


[O///]5007 


15410.8+1.9 


2.0771 +0.0004 


(133.1±7.9)E-17 


72.0E-19 


>90 


(43.0±2.5)E+42 


[JV//]6548 


20157.2+?? 




(5.0±1.0)E-17 


34.9E-19 


>7 


(1.6±0.3)E+42 


Ha 


20203.1+4.9 


2.0776 ± 0.0007 


(58.6±4.1)E-17 


34.9E-19 


>80 


(18.9±1.3)E+42 


[A7/]6583 


20265.9+7.8 




(14.9+3. l)E-17f 


34.9E-19 


>20 


(4.8+1. 0)E+42 


BX1 


Lya(R) 


3921.8+0.6 


2.2260 ± 0.0005 


(10.1±0.8)E-17 


4.9E-19 


150 


(3.9±0.3)E+42 


H/3 


15691.5+28.8 




(1.6±1.5)E-17 


10.3E-19 


>8 


(0.6±0.6)E+42 


[0///J4959 


15986.7+5.7 




(4.3±0.4)E-17 


10.3E-19 


>20 


(1.6±0.2)E+42 


[0///J5OO7 


16149.8+1.9 


2.2246 ± 0.0004 


(12.7±1.2)E-17 


10.3E-19 


>60 


(4.8±0.5)E+42 


[A7/]6548 


21110.3+47.8 




(0.5±1.0)E-17 


7.0E-19 


>3 


(0.2±0.3)E+42 


Ha 


21158.3+5.4 


2.2231 + 0.0008 


(4.1±0.9)E-17 


7.0E-19 


>28 


(1.6±0.3)E+42 


[A7/]6583 


21224.1+?? 




(1.6±2.8)E-17 / 


7.0E-19 


>11 


(0.6±2.8)E+42 


BX10 


Lya(R) 


4403.4+0.3 


2.622 ± 0.0002 


(9.7±0.5)E-17 


7.3E-19 


103 


(5.5±0.3)E+42 


H/3 


17575.0+19.9 




(7.0±6.3)E-17 


44.6E-19 


>8 


(4.0±0.4)E+42 


[0///]4959 


17922.0+20.6 




(6.5±2.4)E-17 


44.6E-19 


>7 


(3.7±1.4)E+42 


[O///]5007 


18095.5+11.6 


2.6132 + 0.0026 


(19.3±7.0)E-17' 


44.6E-19 


>22 


(11.0±4.0)E+42 


Ha 


6562.82x3.614 




<17.5E-17 


35.3E-19 




<9.9E+42 


BX4 


[0///J5OO7 


5007.0x3.505 




(<10.0)E-17' 


18.9E-19 




<5.0E+42 


[JV//]6548 


22948.4+?? 




(1.9±1.0)E-17 


20. IE- 19 


>5 


(0.9±0.5)E+42 


Ha 


23000.9+8.7 


2.5038 ± 0.0013 


(16.6±3.3)E-17 


20. IE- 19 


>40 


(8.2±1.6)E+42 


[A7/]6583 


23071.0+10.6 




(5.8±3.0)E-17 


20.1E-19 


>14 


(2.9±1.5)E+42 


BX7 


H/3 


4863.0x3.418 




<7.5E-17 


15.4E-19 




<3.4E+42 


[0///]4959 


16949.7+10.0 




(4.4±0.7)E-17 


15.4E-19 


>13 


(2.0±0.3)E+42 


[O///]5007 


17113.8+3.6 


2.4171 + 0.0007 


(13.0±2.4)E-17 / 


15.4E-19 


>40 


(6.0±1.0)E+42 


[A7/]6548 


22362.9+82.1 




(0.5±0.5)E-17 


10.7E-19 


>1 


(0.2±0.2)E+42 


Ha 


22413.8+4.0 


2.4143 + 0.0006 


(10.5±1.6)E-17 


10.7E-19 


>50 


(4.8±0.7)E+42 


[A7/]6583 


22483.4+32.7 




(1.6±1.6)E-17 / 


10.7E-19 


>7 


(0.7±0.7)E+42 



Notes. J means sky line residuals can affect the flux estimation. (R) refers to the fact that the rest-frame UV spectrum was normalized to the galaxy 
R band magnitude. Lya emission line fluxes of 'BX' galaxies are estimated from GOODS public spectra (5.7 A/pixel, VIMOS-LR); of LAE27 are 
calculated from spectra obtained with VIMOS (2.57 A/pixel, VIMOS-MR), of z3LAE2 and LBG1 1 with IMACS (1.3 A/pixel) by MUSYC. Optical 
emission lines are obtained from our 7 A/pixel spectra. To determine redshifts we assumed vacuum X Lya = 1215.67 A, A[oni]5!m = 5008.239 A, 
J-h« = 6564.614 A. Column 1, 2, and 3: emission line names, their central wavelengths, and implied redshift. An '??' indicates a formally infinite 
error. Column 4: emission line integrated fluxes. Column 5: continuum flux densities. Column 6: obs-frame equivalent widths of the emission lintl 
listed in the 1st column. In the case of lower limits, the equivalent width is estimated from the 2tr upper limit of the continuum. The continuum 
<t is defined as its root-mean-square. Column 7: emission line luminosities. For a few of the DB sources we report photometric Ha luminosity as 

phot 
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Table 4: Continued: Emission line fluxes for all the sources 



Line 


^obs- frame 
central 


Zspec 


Flux 


rms ix (continuum) 


EW T obs- frame 


L 




A 




erg sec" 1 cm" 2 


erg sec"' cm" 2 A" 1 


A 


erg sec" 1 


BX14 


Up 


16373.0±9.6 




(3.3±2.6)E-17 / 


20.8E-19 


>15 


(1.4±l.l)E+42 


[0///]4959 


16696.2+11.7 




(5.2±1.0)E-17^ 


20.8E-19 


>24 


(2.2±0.4)E+42 


[O///]5007 


16857.8+3.9 


2.3660 ± 0.0008 


(15.6±3.0)E-17 


20.8E-19 


>70 


(6.7±1.3)E+42 


[JV//]6548 
Ho- 


22049.7±?? 
22099.9±4.1 


2.3665 ± 0.0006 


(0.2±0.6)E-17' 
(13.3±2.3)E-17 


14.0E-19 
14.0E-19 


>0.7 
>47 


(0.08±0.3)E+42 
(5.7±1.0)E+42 


[N//]6583 


22168.6±51.6 




(0.6±1.9)E-17^ 


14.0E-19 


>2 


(0.3±0.8)E+42 


BX16 


Lyo-(R) 


4019.0±1.2 


2.3060 + 0.0010 


-(6.3+1. 1)E-17 


12.3E-19 


-45 




up 


16129.1+?? 




(2.9±1.3)E-17 


13.0E-19 


>20 


(2.4±0.7)E+42 


[0///]4959 
[O///]5007 


16447.5+15.0 
16606.7+5.7 


2.3159 + 0.0011 


(1.9±0.6)E-17^ 
(5.6±1.7)E-17 


13.0E-19 
13.0E-19 


>6 
>40 


(0.8±0.2)E+42 
(1.2±0.5)E+42 


[JV//]6548 


21718.1+?? 




(0.6±0.4)E-17^ 


11.5E-19 


>3 


(0.3±0.2)E+42 


Ho- 


21767.5+4.3 


2.3159 + 0.0007 


(10.4±1.6)E-17 


11.5E-19 


>45 


(4.3±0.7)E+42 


[N//]6583 


21835.2+13.4 




(\.l±l.A)E-\lf 


11.5E-19 


>7 


(0.7±0.6)E+42 


LBG11 


Lya(R) 


5112.4+0.6 


3.2054 ± 0.0005 


(3.0±0.7)E-17 


11.6E-19 


34 


(2.8±0.7)E+42 


W 
[0///]4959 
[O///]5007 


20432.6+14.4 
20836.0+16.0 
21037.7+5.0 


3.2006 + 0.0010 


(2.9±2.4)E-17 
(2.6±0.8)E-17' 
(8.0±2.3)E-17 


8.2E-19 
8.2E-19 
8.2E-19 


>18 
>16 
>50 


(2.7±2.2)E+42 
(2.4±0.7)E+42 
(7.5±2.1)E+42 


LAE27 


Lyo-(R) 


4966.7+0.3 


3.0856 ± 0.0002 


(2.2±0.2)E-17 


4.4E-19 


105 


(1.9±0.2)E+42 


up 

[0///]4959 
[O///]5007 


19860.8+48.8 
20252.8+16.4 
20448.9+5.0 


3.0830 ± 0.0001 


(12.0±13.3)E-17 
(8.2±2.6)E-17^ 
(25.8±7.7)E-17 


32.5E-19 
32.5E-19 
32.5E-19 


>19 
>14 
>41 


(10.1±11.3)E+42 
(7.0±2.2)E+42 
(21.9±6.5)E+42 


z3LAE2 


Lya(R) 


5002.4+0.5 


3.1149 + 0.0004 


(2.3±0.5)E-17 


9.9E-19 


74 


(2.0±0.4)E+42 


Up 


20000.5+11.4 




(12.3±5.0)E-17' 


29.6E-19 


>20 


(10.7±6.8)E+42 


[0///]4959 
[O///]5007 


20395.3+11.1 
20592.7+3.1 


3.1118 + 0.0006 


(7.0±1.5)E-17^ 
(20.9±4.4)E-17 


29.6E-19 
29.6E-19 


>12 
>35 


(6.1±1.3)E+42 
(18.2±3.8)E+42 


z2LAE2 


Up 


14849.6+5.7 




(11.7±4.9)E-17^ 


20.0E-19 


>30 


(3.7±1.6)E+42 


[0///]4959 
[O///]5007 


15142.8+10.3 
15289.3+4.4 


2.0528 ± 0.0009 


(11.4±1.5)E-17 
(34.1±4.6)E-17 


20.0E-19 
20.0E-19 


>30 
>90 


(3.6±0.5)E+42 
(10.8±1.5)E+42 


Ho- 


6562.82x3.054 




<15.0E-17' 


28.0E-19 




<4.8E+42 


DB8 


Lya 












<2.8E+41'""" 


up 


15952.8+30.9 




(4.3±3.4)E-17 


15.0E-19 


>14 


(1.8±1.4)E+42 


[0///]4959 
[O///]5007 


16288.1+28.2 
16445.8+8.0 


2.2837 ± 0.0015 


(4.1±0.9)E-17 / 
(12.5±2.8)E-17 


15.0E-19 
15.0E-19 


>15 
>45 


(1.7±0.4)E+42 
(5.1 + 1. l)E+42 


Ho- 


6562.82x3.285 




<5.0E-17 


9.5E-19 




2.3E+41'""" 


DB12 


Lya 












1.6E+42"' M 


up 


15518.0+10.1 




(4.1±3.0)E-17 


16.0E-19 


>12 


(1.5+1. l)E+42 


[0///]4959 
[O///]5007 


15824.3+29.3 
15977.5+7.0 


2.1902 + 0.0014 


(4.0±1.4)E-17 / 
(12.0±4.0)E-17 


16.0E-19 
16.0E-19 


>12 
>27 


(1.5±0.5)E+42 
(4.4±1.5)E+42 


[N//]6548 
Ha 


20907.8+10.4 
20955.4+8.0 


2.1922 + 0.0012 


(1.3±0.9)E-17 
(18.1±3.2)E-17 


18.3E-19 
18.3E-19 


>3 
>50 


(0.5±0.3)E+42 
(6.7±1.2)E+42 


[N//]6583 


21020.5+33.0 




(3.9±2.6)E-17 


18.3E-19 


>10 


(1.4±0.9)E+42 


DB22 


Lya 












<2.8E+41'""" 


up 


4863.0x3.212 




<12.5E-17 


22.8E-19 




<4.9E+42 


[0///]4959 


15930.3+34.4 




(6.9±1.7)E-17 


22.8E-19 


>14 


(2.7±0.7)E+42 


[O///]5007 


16084.5+11.8 


2.2116 + 0.0024 


(20.7±5.7)E-17 


22.8E-19 


>40 


(8.0±2.2)E+42 


Ho- 


6562.82x3.212 




<6.0E-17' 


13.1E-19 




3.9E+41"' M 
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Appendix A: MMIRS data reduction and 
performance 

A.1. Reduction procedure 

The raw mask frames, output of MMIRS, are multi-extension 
files. The 3nd to Nth extensions are characterized by individ- 
ual exposure times equal to t_ramp up to t_total-t_ramp, where 
t_ramp is the ramp-up time and t_total is the total exposure time 
of the frame. The read-out is done using a ramp-up procedure, 
optimized to avoid saturation. It is made every t_ramp seconds, 
t_total/t_ramp times. The 1st extension in the multi-extension 
MMIRS frame contains the total exposure time, while the 2nd 
contains the zero-level counts. 

Dark frames are obtained in the morning for every kind 
of exposure taken during the previous night and in series of 
5 (ramped-up multi-extension darks should be subtracted from 
the ramped-up science frames). We average combined multi- 
extension dark frames, using IRAF msc red. combine task. 

The following step of our reduction pipeline is collapsing all 
of the information enclosed in the multi-extension files. We used 
mmfixen code. The information contained in the N extensions 
are combined, saturated pixels are filled and cosmic rays are re- 
moved. The output is a new multi-extension file, in which the 1st 
extension is the science frame normalized to counts/sec, the 2nd, 
3rd and 4th contain statistical information about the image and 
the fitting process. 

We ran mmfixen code for all the ramped-up science frames 
and then we extracted their 1st extensions. In the case of 
t_total<t_ramp we just subtracted the zero-level counts (2nd ex- 
tension) to the 1st extension. It happened for the majority of flat 
frames and for standard stars. We then normalized flat frames, 
slit by slit, using IRAF twodspec.longslit. response task and we 
divided science and standard frames by the normalized flat. 

COSMOS software was used to generate the wavelength 
map needed for wavelength calibration. We first obtained a 
rough wavelength map using the calibration lamp. The spectra 
of an Argon 9 lamp are obtained through each slit of a mask. 
Lamp emission line pixels were related to known wavelengths 
using COSMOS align-mask and adjust-off set codes. An im- 
proved wavelength solution and spectral map were calculated 
using a few isolated sky lines in science frames. 10 We then ran 
COSMOS ad just-map to obtain the final solution, which con- 
tains the best wavelength map for all the mask slits together and 
it is satisfied by the entire series of sky lines from ~ 14000 to 
-22000 A. 

We applied an ABBA procedure to get rid of prominent sky 
lines. In our May, 2011 run, sky line variations and along indi- 
vidual lines were stronger. Therefore, we got a better sky sub- 
traction in Nov frames. We had the case of the first night of May 
run in which SDSS mask frames were not ramped-up. We just 
dark-subtracted and flat-divided them and directly applied this 
ABBA procedure. We then removed cosmic rays using IRAF 
xdimsum.xzap task to ABBA frames to eliminate cosmic rays, 
coming from A, and to the inverse of ABBA to eliminate cosmic 
rays, coming from B. We checked the alignment of the ABBA 
frames, looking at the continuum spectra of alignment stars (also 
slit stars in the case of May run) and combined all the aligned 
ABBA frames of the same mask. 

We finally applied COSMOS extract-ldspec code to wave- 
length calibrate the completely reduced frame and divide it in 
individual-slit 2D spectra. 

5 MMIRS/MMIRS+Pipeline/Wavelength+Calibration/Argon+Lines 
10 MMIRS/MMIRS+Pipeline/Wavelength+Calibration/OH+Sky+Lines 



A. 2. Flux calibration 

We used IRAF twodspec.apextract.apall task to extract ID 
spectra. For standard stars we located the extraction aperture in 
a region of their continuum not affected by atmospheric absorp- 
tion. For science spectra we chose the pixel position of the few 
prominent emission lines. 

The three standard stars chosen at the observatory are tel- 
luric stars. Therefore they allowed to perform the correction for 
telluric features of science frames and also their flux calibration. 
Spectra coming from the same star, but obtained through dif- 
ferent slits along a mask, cover different wavelength regions, de- 
pending on the position of the slit along the dispersion axis. They 
present comparable shapes in the common wavelength regions, 
but their relative intensity depends on the mask orientation and 
on the sky conditions, that can change with time. These varia- 
tions are reflected in the science frames as well. So we average- 
combined all the ID spectra of the same telluric stars of the run, 
taken at the same airmass of the science frames. 

We extracted ID spectra of the objects in the center of the 
slits as well as the ones of the background either on the upper 
or lower side of the slits to estimate background fluctuations, 
after sky emission line removal. In Fig. A.l we show z2LAE3 
raw(smoothed) spectrum in black(red) and the background spec- 
trum in green. This background spectrum is a proxy of the noise 
spectrum of this object. 

Appendix B: MMIRS 2D spectra and performance 

As we could see in Table 2, z < 1.4 'low-z' galaxies and slit 
stars were located in MMIRS masks as a test of its performance. 
Fig. B.l shows 2D spectra of two slit stars in SDSS (upper) 
and EHDF-S (lower) masks. A star of K AB = 16.42 starts to 
show saturation in an MMIRS spectrum, while the continuum 
of a Kab = 17.81 star is clearly detected. Regions of strong 
atmospheric absorption are clearly seen, at ~ 14000 A (in the 
left side of the spectra), at ~ 19000 A (in the center-right side of 
the spactra). In these specific cases the slits were located at the 
same position along the dispersion axis of MMIRS mask. In Fig. 

B. 2 we present the 2D spectra of 5 low-z galaxies obtained with 
SDSS mask. The top panel shows the source with the brightest 
Kab band mag of 17.65. The continuum appears well detected 
for sources with Kab ^ 19.39. 
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1.2xl0 4 1.4x10* 1.6xl0 4 l.BxlO 4 2xl0 4 2.2xl0 4 2.4x10" 

A[A] 

Fig. A. 1 : ID spectrum of z2LAE3 along the entire MMIRS coverage. The black (red) curve represents the raw (smoothed) spectrum. Grey vertical 
lines indicate the wavelength of the dominant sky emission lines. Cyan curve represents the atmospheric transmission at Las Campanas observatory. 
At A ~ 14000 and ~ 19000 A the transmission is very low. The green curve represents the background spectrum extracted far from the center of 
the slit where the source is positioned. Sky line residuals produce higher uncertainty. Telluric correction fixes the continuum shape but increases 
the uncertainty. Wavelength regions of high noise are not shown in the smoothed spectrum. 



K=17.81 




Fig. B.l: 2D spectra of two slit stars in the SDSS and EHDF-S masks with K AB = 17.81 (up) and K AB = 16.42 (low), respectively. 




Fig. B.2: 2D spectra of the 'low-z', low-redshift galaxies obtained with SDSS mask. The spectra are aligned in wavelength. The top panel shows 
SDSS28 with the highest K AB band mag of 17.65. These sources spanned K AB band magnitudes from 17.65 (SDSS28), 19.39 (SDSS12) to 20.46. 
In the center of the stripes (~ 19000 A), the decrease in flux indicates the wavelength region of strong atmospheric absorption. We can also see 
regions of high Poisson noise due to sky emission lines. 
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